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Abstract 

I briefly review the arguments why the braneworld models with infinite- 
volume extra dimensions could solve the cosmological constant problem, evad- 
ing Weinberg's no-go theorem. Then I discuss in detail the established prop- 
erties of these models, as well as the features which should be studied further 
in order to conclude whether these models can truly solve the problem. This 
article is dedicated to the memory of Ian Kogan. 
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1 Two puzzles: cosmological constant and cosmic 
coincidence 



Perhaps, the least understood problem of particle physics, gravity and cosmology 
is that of the Cosmological Constant. The problem stems from a huge mismatch 
between observational data and theoretical expectations. The problem can be briefly 
outlined as follows. The structure in the Universe (galaxies, clusters etc.) could have 
formed only if the acceleration rate of the expansion i^A is less than the number 
that roughly equals to the present-day value of the Hubble parameter Hq P, 

Hjy <Ho ^ 10-33 eV . (1) 

Recent observations |2] appear to confirm a nonvanishing value of i^A which nearly 
saturates the upper bound in Eq. (0). On the other hand, in general relativity (GR) 
H\ determines the scalar curvature of space-time and is related to the vacuum 
energy density S as follows: 

Mi.Hl^S. (2) 

Here Mpi ~ 10^^ GeV is the Planck mass. Moreover, a natural value of S due to 
zero-point oscillation energies of known elementary particles can be estimated as 
S > (TeV)'^. Substituting this value of £ into (j2)), one finds i^A ^ IQ-^eV, which 
is grossly inconsistent with (0). This is the essence of the cosmological constant 
problem (CCF)} 

Historically, the CCP was formulated long before the discovery of the cosmic 
acceleration. One of the first published works on the subject was by Zel'dovich 
[3] in 1967 where he estimated the contribution of zero-point oscillation energies of 
nucleons to the vacuum energy density. Naturally, he found 

S ~ ~ (GeV)^ 

which already gives a result grossly inconsistent with (P).^ 

The problem only worsened in the 1970s and '80s as particle physics made huge 
steps forward in understanding Nature at exceedingly shorter distances. This only 
increased the value of the maximal momentum accessible to a particle whose zero- 
point oscillation energy contributes to the energy density of the Universe. Therefore, 
estimate of a natural value of £ grew up to £^ > (TeV)"^. Since then, theorists 
continuously worked on the problem. Although no satisfactory solution has been 
found, nevertheless, as it usually happens, many new useful theoretical aspects got 
uncovered in the search process (for a review see, e.g., Q). 

^We only discuss positive £, all the essential arguments apply to negative £ too. 

^According to colleagues who witnessed those developments Zel'dovich got so frustrated by this 
problem that he practically stopped doing particle physics and turned his attention to astrophysics 
where he has made great contributions in the 1970s. I thank Misha Shifman for his recollections 
of that period. 
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A dramatic reshaping of the subject took place with the discovery of the cosmic 
acceleration at the end of the previous millennium [2]. This discovery triggered a 
fresh tremendous interest in the problem and motivated recent developments. So 
far the discovery only sharpened the status of the problem, making us to realize 
that there are two puzzles that we have to face. In a conventional formulation these 
puzzles can be spelled out as follows: 

(i) Why is the vacuum energy in the Universe so much smaller than any reason- 
able estimate that follows from particle theories? This is the "old" CCP. 

(ii) Why is the vacuum energy in the Universe comparable to matter energy? Or, 
do we live in a special epoch when the magnitudes of the above quantities roughly 
coincided This is the so called cosmic coincidence problem (CoCoP). 

A priori one could choose to adjust by hand the renormalized values of the 
vacuum energy density to be equal to £^ ~ (10~^ eV)"^, to make it consistent with 
observations. There are many classical as well as quantum-mechanical contributions 
to the vacuum energy such that (a) some of these contributions are many orders of 
magnitude larger than S ~ (10~^ eV)^ and (b) some of these contributions differ 
from each other by many orders of magnitude. Hence, this adjustment requires an 
incredible fine-tuning of the parameters. Adopting the fine-tuning, we could suc- 
cessfully parametrize the observed cosmological evolution of the Universe. However, 
the fundamental questions (i) and (ii) would still remain open since it is not clear 
why such different contributions to the vacuum energy had to cancel to such an ex- 
traordinary accuracy and why the result of that extraordinary cancellation should 
be of the same order as the present-day value of the matter density in the Universe. 

In spite of numerous attempts, neither of the above puzzles have satisfactory 
explanations so far.^ Most of the approaches that have been developed to solve (i) 
are disfavored by a general no-go theorem formulated by S. Weinberg jl]. As to the 
solution of (ii), it seems more reasonable to think about it only in the context of (i). 

Let us point out that the formulation of the question (i) itself contains a loophole 
which might be suggestive of a new approach to the solution of CCP. Indeed, we 
have no direct experimental way to measure S. Instead, we measure space-time 
curvature through cosmological observations, and then determine S through the 
Einstein equations. Thus, claiming that S should be small we implicitly assume that 
the Einstein equations are valid for arbitrarily large length scales. This assumption 
may or may not be correct. This suggest an alternative approach where S keeps 
its natural value S > TeV^, but laws of gravity are modified so that large vacuum 
energy density does not give rise to large space-time curvature. Since the discrepancy 
between the theory and experiment manifests itself at enormous distances ~ Hq^ ~ 

^We are interested in an explanation in terms of a low-energy theory. The anthropic approach 
(for a review see, e.g., 0) seems to give answers to both (i) and (ii), and is certainly a logical 
possibility. However, this is an orthogonal approach. Another logical possibility is that the problem 
can never be understood in terms of low-energy dynamics and is only solved due to very contrived 
effects of ultraviolet physics, which, in fact, might not be as contrived as it might seem, because 
of symmetries of string theory 
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10^*^ cm (i.e., at extremely low energies), to address CCP it is natural to modify 
gravity in the infrared domain (IR). 

Construction of such models was motivated by the advent of the braneworld 
paradigm, where the standard-model fields are localized on a brane while grav- 
ity propagates in the bulk [7] (for earlier models see jHJ 13 UHl; recent reviews 
can be found in Refs. fT] ^1 El)- However, making a consistent theory of the 
IR-modified gravity became possible only in models with infinite-volume extra di- 
mensions ^1 ! where gravity on the brane transforms from four-dimensional to 
higher- dimensional at very large distances. Historically the first was a proposal of 
Ref. [inj, which gives a brane- world realization of a massless and massive gravity. 
This was followed by an early proposal of a theory of a metastable graviton p!7] . 
However, the latter turned out to be an internally inconsistent theory ^Hl E] • 

The model of Ref. Jl], and its higher-dimensional generalizations ^Hj, paved the 
way to new possibilities of addressing the cosmological constant problem through IR 
modification of gravity, where the vacuum energy (the brane tension) mostly curves 
the bulk, while ordinary gravity is trapped on the brane at observable distances by 
the presence of a large Einstein-Hilbert action localized on the brane. 

A specific proposal along these hues was worked out in Ref. [201 , where it is 
argued that the graviton propagator is modified in the infrared in such a way that 
large wavelength sources, such as the vacuum energy, gravitate very weakly. As a 
result, even a huge vacuum energy does not curve our space. On the other hand, 
short wavelength sources, such as planets, stars, galaxies and clusters gravitate 
(almost) normally. The four- dimensional (4D) nonlocal counterpart with similar 
properties was proposed in Ref. pT| . 

We will discuss the framework of Refs. [T3] where gravity in general, and 
the Friedmann equation, in particular, are modified for wavelengths larger than 
a certain critical value. This setup can evade the Weinberg no-go theorem. The 
cosmological constant problem could then be remedied in the following way: Due to 
the large- distance modification of gravity the energy density £^ > (1 TeV)^ does not 
curve the space as it would do in the conventional Einstein gravity. Therefore, the 
observed space-time curvature is small, despite the fact that £ is huge (as it comes 
out naturally). This is the most crucial point of the approach of Refs. [UIIS] - the 
point where we depart from the previous investigations. Although, as we will see, 
it is still premature to say whether this approach leads to a final solution of CCP, 
nevertheless, it seems that all necessary ingredients are present in the model. Future 
detailed calculations will show whether or not this development is successful. 

Before delving in the issue we would like to mention that the idea of solving the 
cosmological constant problem in theories with extra dimensions and branes has a 
long history (for earliest works see, e.g., [221I2S1)- However, because of Weinberg's 
theorem, the solution is only possible in theories where the extra dimensions have 
infinite (or practically infinite) volume. Why is this so? A brief answer will be 
presented below (more complete discussions are given in Ref. [20] )• 

Recall that if there is a vacuum energy density £ > TeV*^ in a conventional 4D 
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theory then it unavoidably gives rise to the scalar curvature R ~ H\ determined 
by (j2I)- The vacuum energy density ^ is a source of gravity, and, as such, it has to 
curve the space; the only space in 4D theories is the space in which we live. Hence, 
our space is curved according to (j2I), and this is inconsistent with data. However, 
if there are more than four dimensions, £ could curve extra dimensions instead of 
curving our 4D space [221 12^] ■ Consider the following (4 + A^)-dimensional interval: 

ds" = A\y) g^,{x) dx^dx" - B\y) dy^ - C\y) y" dn%,, , (3) 

where /i, z/ = 0, 1, 2, 3, are the indices denoting our 4D world, while y = ^yf + ... + y%, 
and ynS denote extra coordinates. What we measure in our 4D world is the curvature 
invariants of the metric g^y{x). There can exist solutions to the (4 + A^) -dimensional 
Einstein equations in the form of Q where £ affects strongly the extra space, i.e., 
the functions A, B and C, while leaving our 4D space almost intact, with the 4D 
metric g^v{x) remaining almost flat. 

In this case the energy density £ "is spent" totally on curving up the extra space 
rather than on curving our 4D space. The simplest example of this type is a 3-brane 
in six-dimensional space (a local cosmic string) in which case the tension of the 
brane is spent on creating a deficit angle in the bulk, while the brane world-volume 
remains fiat (for a discussion see [21] )• 

Such a brane could be a good place for our 4D world to live. If one could only 
obtain the laws of 4D gravity on a brane in this setup, this would be considered as 
a solution of CCP that takes into account all classical and quantum contributions 
to the cosmological constant! 

The same arguments would apply to higher codimensions. Therefore, the paramount 
goal is to find a mechanism that would enable one to obtain 4D gravity on the brane 
embedded in infinite- volume bulk.^ 

The remainder of this article describes a method of obtaining 4D gravity on a 
brane in infinite- volume extra space. First, in Sect. 2 we formulate a basic model 
[m [T3]. Then we discuss how this model evades Weinberg's theorem. In Sect. 3 
we consider in detail this model in five dimensions. Although the 5D model is 
known a priori to be unfit to solve CCP, nevertheless, it is instructive to study this 
situation in detail. Most of the intricate properties of the 5D model are understood, 
and one can say that the model with appropriate boundary conditions represents 
a consistent theory of a large-distance modification of gravity. In Sect. 4 we turn 
to similar models in more than five dimensions. Here the situation is different. We 
discuss what is known so far about these models and what needs to be done in 
order to conclude whether this approach can solve CCP. Section 5 contains a brief 
summary. 

■*Any compactification of the above setup with the compactification radii smaller than H^^ 
would give rise to a theory of gravity that flows to the conventional GR in the IR. The latter would 
necessarily face Weinberg's no-go theorem, for details see |2()j . 
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2 The origin of the model 



In this section we will formulate the model which was introduced in 5D space-time in 
Ref. and later generalized to > 6 in We closely follow the presentation 
of Ref. 1201 • 

Consider a brane-world model in a space with (asymptotically) flat infinite- 
volume N extra dimensions. Assume that all known standard-model (SM) particles 
are localized on the brane and obey the conventional 4D laws of gauge interactions up 
to very high energies, of the order of the GUT scale, for instance. The gravitational 
sector, on the other hand, is spread over the whole (4 + A^)-dimensional space. The 
low-energy action of the model is written as 

S = Ml^"^ j d'xd'^y^gn.+Mig) 

+ j d^x^ (£ + Ml^ R + £sm(^, Msm)) • (4) 

Let us discuss various parts and parameters of the action (jH). £sm is the Lagrangian 
for particle physics including all SM fields The parameter Msm denotes the 
ultraviolet (UV) cutoff of SM. Up to that scale SM obeys the conventional 4D laws. 
In the present approach Msm ^ TeV ^ M^. Moreover, cjab stands for a (4 + A^)- 
dimensional graviton {A, B = 0, 1, 2, 3 + A^), while ?/„, n = 4, 5, .., 4 + A^, denote 
"perpendicular" to the brane coordinates. For simplicity we do not consider brane 
fluctuations ^. Thus, the induced metric on the brane is given by 

9,^uix) = g^,uix,yn = 0) . (5) 

Since we discard the brane fluctuations, the brane can be thought, in the 5D case, as 
a boundary of the extra space or an orbifold fixed point (in that case the Gibbons- 
Hawking surface term is implied in the action hereafter). The brane tension is 
denoted by £. 

The first term in is the bulk Einstein-Hilbert action for (4 + A^)-dimensional 
gravity, with the fundamental scale M*. The expression in has to be under- 
stood as an effective low-energy action valid for graviton momenta smaller than 
M^,. Therefore, in what follows we will imply the presence of an infinite number of 
gauge-invariant high-dimensional bulk operators suppressed by powers of M^. 

The second term in (j^ describes the 4D Einstein-Hilbert (EH) term of the 
induced metric. This term plays the crucial role. It ensures that at observable 
distances on the brane the laws of 4D gravity are reproduced in spite of the fact that 

^For notational simplicity we use the convention that the particles physics theory, including 
any grand unification (GUT), possibly SUSY GUT, or any other extension of standard model, is 
denoted as SM. 

^This limitation could be readily lifted. Indeed, including the brane fluctuations would produce 
an almost sterile Nambu-Goldstone boson, and heavy modes which could manifest themselves only 
through generation of an extrinsic curvature term on the brane. 
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there is no localized zero-mode graviton. Its coefficient Mind is another parameter of 
the model. Thus, the low-energy action as it stands is governed by three parameters 
M^,, Mind and E. Let us discuss their natural values separately. 

The parameter Mind gets induced by SM-particle loops localized on the brane. 
Such corrections are cut-off by the rigidity scale of SM, Msm, i-e., the scale above 
which the SM propagators become soft. In the present approach this scale is taken 
to be very high, ^ TeV. In particular, we will take this scale to be comparable with 
the GUT or 4D Planck scale. ^ The loops induce^ the Einstein-Hilbert term in (j^, 

MLv^^(^?), (6) 
where the value of the induced constant Mind is determined by the relation |2Zl I2H] , 

Mfnd = ijd^x {T{x) T(0))/96 . 

The parameter Mind is proportional to the scale Msm and to the number of the SM 
particles.^ Since there are about 60 particles in the Weinberg-Salam model, and 
more are expected in GUT's, the value of Mind should be somewhat larger than 
Msm- In fact, below we define the 4D Planck mass as being completely determined 
by Mind, 

Mpi = Mind ■ (7) 

Thus, the Planck mass is not a fundamental constant in our approach but rather 
a derived scale. We see that the SM loop corrections are capable of creating the 
hierarchy Mind/M*, even if the initial value of Mind/M,,, was not that large. This 
hierarchy does not amount to fine tuning, since such a separation of scales is stable 
under quantum corrections. Indeed, say, M^, gets renormalized by all possible bulk 
quantum gravity loops. However, there are no SM particles in the bulk the only 
scale in there is M^. Therefore, any bulk loop gets cut-off at the scale M*, as it is 
the fundamental gravity scale. While, as we discussed above, the brane SM loops 
are cut-off by the higher scale Msm, and this gives rise to the huge value of Mind on 
the brane. 

Finally, let us discuss the value of the brane tension (brane cosmological con- 
stant), and of the bulk cosmological constant. To this end, we have to specify our 
assumptions regarding supersjTiimetry. We assume that the high-dimensional the- 
ory is supersymmetric, and that supersymmetry is spontaneously broken only on the 

^We set the thickness of the brane A to be determined by the SM scale, A ^ -^sm • This 
might seem a bit unnatural at a first sight, but there are field theory [221 as well as string theory 
constructions [23 of branes where such a "dynamical" width is possible. 

® A/ind can certainly contain as well the tree-level terms if these are present in the original action 
in the first place. We will not discriminate between these and induced terms. Mind will be regarded 
as a parameter that stands in 

^The scalars and fermions contribute to Mjnd with positive sign while the gauge fields with 
negative sign. 
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Figure 1: (A) The world- volume vacuum diagram of the SM fields that renormalizes 
the brane tension (i.e., the 4D cosmological constant). These contributions are 
protected by A/" = 1 world- volume supersymmetry. Therefore, they are cutoff by the 
world- volume SUSY breaking scale Mgusy (B) The world- volume two-point diagram 
that renormalizes (induces) the EH term on the brane. These contributions are 
not protected by A/" = 1 supersymmetry. They can only be protected by conformal 
invariance which in our model is broken at the scale Msu that is close to Mpi. Hence, 
there can be a hierarchy between (B) and (A). (C) The bulk vacuum diagram. Only 
the bulk particles (which do not include the SM particles) are running in this loop 
(bulk particles are denoted by double lines). This diagram is protected by unbroken 
bulk SUSY. Therefore, the cosmological constant in the bulk is zero. (D) The bulk 
two-point diagram which renormalizes the bulk EH term. As in (C), only the bulk 
particles are running in the loop. This diagram is cutoff by the bulk scale M*. 
Therefore, the natural value of the constant in front of the bulk EH term is M^"*"^; 
there is huge hierarchy between this coefficient and that of the world-volume EH 
term coming from (B). 

brane (such a scenario with a non-BPS brane was considered in fTU^). The absence 
of breaking of supersymmetry in the bulk is only possible due to infinite volume 
of the extra space; SUSY breaking is not transmitted from the brane into the bulk 
since the breaking effects are suppressed by an infinite volume factor. Then, the 
bulk cosmological term can be set to zero, without any fine-tuning. On the other 
hand, the natural value of S can be as low as TeV^, since the brane tension can be 
protected above this value by A/" = 1 supersymmetry (note that Mind can only be 
protected by a conformal invariance which we assume is broken at the scale Msm)- 
All these properties are summarized in Fig. 1. 

Let us now turn to the gravitational dynamics on the brane. This dynamics is 
quite peculiar. Despite the fact that the volume of extra space is infinite, an observer 
on the brane measures 4D gravitational interaction up to some large cosmologically 

^°In general, local SUSY in the bulk does not preclude a negative vacuum energy density of the 
order of M^"*"^. However, the latter can be forbidden by an unbroken R symmetry in the bulk. 
Such a symmetry is often provided by string theory. 
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scales. The fact that this is so will be studied in detail in the following sections. 

The no-go arguments discussed in the previous section are not applicable to the 
theories with infinite-volume extra dimensions. The crucial property of this class of 
theories is that despite the unbroken 4D general covariance, there is no 4D zero-mode 
graviton. 4D gravity on the brane is mediated by a collective mode which cannot 
be reduced to any 4D state. The fact most important for us is that the 4D general 
covariance does not require now all states to couple universally to our "graviton" . 
As a result, there is no universal agent that could mediate supersymmetry breaking 
from SM to all existing states. Such a situation is impossible in the finite-volume 
theories where 4D gravity is mediated by a normalizable zero mode, which, by 
general covariance, must couple universally and, hence, mediates supersymmetry 
breaking. Moreover, the effect of the brane cosmological term is to curve the extra 
space without inducing a large 4D curvature. We stress again that this is impossible 
in finite- volume theories (i.e., the theories in which the size of the extra space is 
smaller than the Hubble size -f^o"^) because there the extra components of the metric 
are always heavier than Hq. 

The effective field theory arguments are based on the assumption that there is a 
finite number of 4D degrees of freedom below the scale of the cosmological constant 
that one wants to neutralize. This condition is not satisfied in the present model - 
it is a genuinely high- dimensional theory in the far infrared. Therefore, there is an 
infinite number of degrees of freedom below any nonzero energy scale. As a result, 
there is no scale below which extra dimensions can be integrated out and the theory 
reduced to a local 4D field theory with a finite number of degrees of freedom. In order 
to rewrite the model at hand as a theory of a single 4D graviton, at any given scale, 
we have to integrate out an infinite number of lighter modes. As usually happens in 
field theory, integrating out the light states we get nonlocal interactions. Therefore, 
the resulting model, rewritten as a theory of a 4D graviton, will contain generally- 
covariant but nonlocal terms. The latter dominate the action in the far infrared. Of 
course, in actuality, the full theory is local - the apparent nonlocality is an artifact 
of integrating out light modes. It tells us that a local (4 + A^)-dimensional theory 
can be imitated by a nonlocal 4D model. The nonlocal terms modify the effective 
4D equations and neutralize a large cosmological constant. 

3 Five-dimensional model 

In this section we concentrate our attention on a 5D model of brane induced gravity 
- the so called DGP (Dvali-Gabadadze-Porrati) model jT3]. This model, as a field 
theory, exhibits many unusual and exciting properties that one encounters in theories 
of large distance modified gravity. In a class of Lorentz-invariant theories the model 
is right now the only internally consistent theory of large distance modification of 
gravity (along with some of its higher- dimensional generalizations, see below). 
However, the 5D DGP model cannot solve the "big" cosmological constant prob- 
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lem because the brane in it has only one codimension. In the best case, the 5D 
model can successfully parametrize the accelerated Universe [221 ISHl HI] (although 
the viability of the latter assertion still needs to be established in greater detail, see 
below) . 

Therefore, the 5D DGP model can be regarded as a toy example on which many 
intricate features of large distance modified gravity can be understood. Some of 
these features will be important in searching for the theory of large distance modified 
gravity that could truly solve the "big" cosmological constant problem. 

After reviewing the 5D model in this section, we turn to the more general models 
in D > 6 in the next section. The action of the 5D model jT3] is 

S = Ml, j d^x^R{g) +Ml J d'xdy^n,{g) , (8) 

where R and TZ^ are the four-dimensional and five-dimensional Ricci scalars, respec- 
tively, and stands for the gravitational scale of the bulk theory. The analog of 
the graviton mass is rric = 2M^/M|j. The higher-dimensional and four- dimensional 
metric tensors are related as 

g{x,y = 0) = g{x) . (9) 

There is a boundary (a brane) at y = and Z2 symmetry across the boundary 
is imposed. The presence of the boundary Gibbons-Hawking term is implied to 
warrant the correct Einstein equations in the bulk. Matter fields are assumed to be 
localized on a brane and at low energies, that we observe, they do not escape into 
the bulk. Hence, the matter action is completely four-dimensional Sm = I d'^xLM- 

Our conventions are as follows: tiab = diag [H ] ; A, B = 0,1,2,3, 5 ; fi,^ = 

0,1,2,3. 



3.1 Perturbative expansion in Newton's constant 

A simplest exercise that tells us a lot about the model is to calculate the Green's 
function D^^'""^ and the amplitude of interaction of two sources T^^ on the brane 

•^1— graviton — ^ [iv ^ -'- a(3 ■ V-'-"/ 

In order to perform perturbative calculations one has to fix a gauge. One choice, 
that was adopted in [H], is harmonic gauge in the bulk O^Hab = f^B^S/^- Then, 
the momentum-space one-graviton exchange amplitude on the brane takes the form: 

exp {—p\y\) 

A-graviton(P, y) = 5— , (H) 

+ rricP 

where we denote the Euclidean four-momentum squared by p^, 

2 u 2i2i2i2 2i2i2i2 

P = -P^Pf. = -Po + Pi +P2 +P3 = P4 + Pi +P2 +P3, (12) 
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and 



^ 8vrG;v(r^l - It-T) . (13) 



1/3 = l^V^ ~ 3 

In the expressions above p stands for the square root of p 



The euchdean amphtude (fTT|) was constructed by imposing the decreasing boundary 
conditions in the y direction. 

In principle, one could choose the other sign of the square root while solving the 
equations and obtain the euclidean amplitude that grows with y 

T?/3 exp (p|?/|) 

-4i_graviton(P, V) = • (15) 

p^ — nicP 

The latter expression differs from the one in (jllj) not only in its y dependence, but 
also by the position of the pole in the denominator. 

The above two solutions ()11|) and ()15|) are distinguished from each other by 
the choice of the boundary conditions at y — ±00. The choice of the decreasing 
boundary condition in (lllj) is conventional, and as we will see below, under this 
choice one obtains the expected results - the 4D gravity at r -C Tc = is 
smoothly transitioning to 5D gravity at r ^ Tc. On the other hand, the choice of 
the growing boundary conditions in ()15|) might seem somewhat unusual. However, 
as we will see below, the Minkowski space is unstable for this choice, and as a result 
one obtains the so called selfaccelerated space [22] which can be used to describe 
the accelerated expansion of the Universe fHO j. 

To reveal these properties we study the pole structure of (fTT|) and (fT^. Let us 
start with ()11|). We refer to the branch with this choice of the boundary conditions 
as the "conventional branch" as opposed to the "selfaccelerated branch" specified 
by (|15p. The equation determining the poles on the conventional branch is 

2 



p + rricyp'^ = 0. (16) 

Hence there are at least two poles, one at = and another one at p^ = —rric p. Our 
goal is to establish where this poles are located on the complex plane of minkowskian 
momentum square p^. The transition between euclidean momentum square and 
the minkowskian momentum square is as follows: 

/ = e->J. (17) 
Using this we find poles in minkowskian momentum square 

pI = 0, pl = mle-'-. (18) 
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As it can be checked, the residue of the = pole is zero. Therefore, there 
is no massless mode that can mediate interactions in this model. The remaining 
pole is located on a nonphysical Riemann sheet, pointing to a resonance nature of 
the graviton. The residue in this poles can also be calculated and it is positive - 
corresponding to a residues of a positive norm state. Hence, on the conventional 
branch we obtain one metastable graviton with the lifetime r ~ ~ -f^o"^- 

Let us now turn to the "selfaccelerated branch" (jl5|) . The poles of this expression 
are now determined by 

— mc\ljP' = 0. (19) 

Using the same arguments as above we find the poles, 

pj = 0, pl = my-. (20) 

As before, the pole at = has zero residue, hence there is no massless graviton 
in this case either. On the other hand, the pole at = m^e*'^ has a positive residue 
of a positive norm state. However, this pole is located on a physical Riemann sheet. 
Therefore, it describes a tachyon-like state. This signals that the Minkowski space 
is unstable on this branch. The instability should grow with time as 

6"^=* . (21) 

This is a welcome feature since this instability could signal that the background 
should be readjusted and that the curvature of the new background should be of 
the order of m^. On the other hand, this is roughly the curvature that is needed to 
describe the accelerated universe (for a modern review on theory and observations, 
see (SI). 

The above perturbative arguments can be generalized to a full-fledged nonpertur- 
bative analysis by looking at exact cosmological solutions of the model [211 E3 ISO] • 
One can the exact cosmological equations for studying the evolution on the brane. 
What is important here is the expression for the Friedmann equation on the brane. 
In terms of the Hubble parameter of the 4D brane world- volume H, the latter equa- 
tion takes the form 

(I) Conventional branch, i.e., decreasing boundary conditions at y —>■ ±oo (com- 
pare with /f7^)): 

+ nicH = 0. (22) 

There are two solutions to the above equation: 
Solution (A) 

H = 0. (23) 

This solution corresponds to the Minkowski space of the conventional branch. Small 
perturbations about this space are stable. On this solution the cosmological evolu- 
tion transitions from a 4D regime when H ^ rric to the 5D regime when H <C rric- 
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This behavior might be useful for certain cosmological issues, however, it cannot 
explain the accelerated expansion of the Universe. 
Solution (B) 

H = -rUc . (24) 

This corresponds to a collapsing Universe with the scale factor exp(— met) and the 
typical time scale determined by re- 

(II) Self accelerated branch, i.e., increasing boundary conditions at y ^ ±cxd 
(compare with (03)): 

- m^H = 0. (25) 

These are empty space Friedmann equations. 
Solution A' 

H = 0. (26) 

This is a Minkowski solution of the selfaccelerated branch. However, as we discussed 
above, the small perturbations about this branch reveal the exponential instabilities 
of the type pij) with the typical time scale determined by Tc- There are two questions 
in this regards: 

(1) Where this instability leads the theory? 

(2) Whether this instability can be used to mimic the accelerated expansion of 
the Universe? These questions were not studied yet. 

Solution B' 

H = rUc. (27) 

This is a selfaccelerated solution found by Deffayet [23 • This solution was shown to 
describe successfully the accelerated expansion of the Universe jSUlEIl- The question 
whether this solution itself is stable with respect to small fluctuations needs further 
detailed studies in the light of the results of Refs. [3311221 where it was shown that in 
a particular limit of the theory there is a ghost-type excitation on the selfaccelerated 
background (P7j). The question whether this ghost is present on the selfaccelerated 
background in the full theory and is not an artifact of the particular limit taken in 
Refs. [HH iSSj needs to be studied. 

Some of the discussions presented above were based on purely perturbative ar- 
guments (although in all the cases the exact results could also be obtained). In 
this regard, it is appropriate to wonder about the limitations of the perturbation 
theory in the present case. As we will see, it turns out that the naive perturbative 
expansion in Newton's constant breaks down unusually early as compared to the 
standard GR. In general terms the reason for this breakdown is as follows. The 5D 
model has two dimensionful parameters: the Newton constant Gj^ and the graviton 
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lifetime m^. The naive perturbative expansion in powers of Gn is contaminated by 
powers of l/rric. Hence, for small values of perturbation theory breaks down for 
the unusually low value of the energy scale. 

The reason for the breakdown of perturbation theory at a low scale can be traced 
back to terms in the graviton propagator that contain products of the structure 

, (28) 

rricP 

with similar structures or with the flat space metric. These terms do not manifest 
themselves in physical amplitudes at the linear level since they are multiplied by 
conserved currents, however, they enter nonlinear diagrams leading to the break- 
down of perturbation theory similar to massive non-Abelian gauge fields or massive 
gravity. However, this breakdown is an artifact of an ill-defined perturbative ex- 
pansion - the known exact solutions of the model have no trace of the breakdown 
(see Refs. jSHllini)- This shows that if one sums up all the tree- level perturbative 
diagrams, then the breakdown problem should disappear. 

For a source of mass M and the Schwarzschild radius tm = 2GnM, the pertur- 
bative breakdown scale takes the form [IHI 

= (rMr2)V3. (29) 

This is a scale at which nonlinear interactions in a naive perturbative expansion in 
Gn become comparable with the linear terms (below we will discuss in detail the 
physical meaning of this scale). For a source such as the Sun, the hierarchy of the 
scales is as follows: 

tm ~ 3 10^ cm < r, ~ 3 10^° cm < ~ 10^^ cm . (30) 

It is interesting to note that for cosmological solutions of the Friedmann-Robert son- 
Walker (FRW) type r^, ~ Vc- (The same is true for very law energy density sources.) 
Therefore, the perturbative calculations described above give valid results for the 
FRW type solutions only at distance/time scales larger that r^. This is confirmed 
by exact cosmological solutions. 

Hence, the conclusions of the above perturbative calculations can be used to 
state that there are two branches of solutions, that have different behavior at t > Tc. 
These are the conventional solution and the selfaccelerated branch described above. 
The Minkowski space is stable on the conventional branch, however it is unstable 
on the selfaccelerated branch. 



3.2 Constrained perturbation theory 

In this section we consider a possibility of modifying the linearized perturbation 
theory in the DGP model by introducing certain new terms that would enable to 
remove the singular in rric terms from the propagator. We closely follow Ref. |¥T] . 

^^Unfortunately, the massive gravity in 4D [201 is an unstable theory |37] with an instability 
time scale that can be rather short |38|. 
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For this we recall that the the breakdown of the perturbative expansion in 
can be traced back to the expression for the trace of h^^, = g^^ — r]^^ which in the 
harmonic gauge takes the form [T^ 



T 

SnicP 



h^,{p,y = 0) = -TZ—^- (31) 



(Hereafter the tilde denotes Fourier-transformed quantities, and we put SuGn = 
1.). From this expression we learn that: (i) /i^ is a propagating field in this gauge; 

(ii) h'^ propagates as a 5D field, i.e., it does not see the brane kinetic term; (iii) 
The expression for /i^ is singular in the limit rric 0. The gauge dependent part of 
the momentum-space propagator D{p,y) contains the terms p^p^h, which, due to 
fj31|) . give rise to the singular in l/m^ term. Hence, to understand the origin of the 
breakdown of perturbation theory, one should look at the origin of the l/m^ scaling 
in dSD). 

The singular behavior of is a direct consequence of the fact that the four- 
dimensional Ricci curvature R{g) in the linearized approximation is forced to be 
zero by the {55} and/or {[xb} equations of motion. This can be seen by direct 
calculation of R and of those equations, but it is more instructive to see this by 
using the Arnowitt-Deser-Misner (ADM) decomposition. The {55} equation reads 

R = {K':,f - Kl, , (32) 

where K^y denotes the extrinsic curvature. Since K ~ 0{h) the above equation 
implies that the four-dimensional curvature R ~ 0{h'^) and in the linearized order 
R vanishes. Let us now see how this leads to the singular behavior of h in (jHlj) . 
The junction condition across the brane contains two types of terms: there are 
terms proportional to ttIc and there are terms that are independent of nic. The 
former come from the bulk Einstein — Hilbert action while the latter appear due 
to the world-volume Einstein-Hilbert term. In the trace of the junction condition 
the TTic independent term is proportional to the four- dimensional Ricci scalar R. 
On the other hand, as we argued above, R has no linear in h term in the weak- 
field expansion, simply because these terms cancel out due to the {55} and/or 
{/i5} equations. Therefore, in the linearized approximation the junction condition 
contains only the terms that come from the bulk. These terms are proportional to 
mji. This inevitably leads to the trace of h (j31|) that is singular in the uic 
limit and triggers the breakdown of the perturbative approach as discussed above. 

The above arguments suggest that the two limiting procedures, first truncating 
the small h expansion and only then taking the rric — > limit, do not commute 
with each other. Therefore, the right way to perform the calculations is either to 
look at exact solutions of classical equations of motion, as was argued in jSHllini, or 
to retain at least quadratic terms in the equations. The obtained results won't be 
singular in the rric ~^ limit. 
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However, neither of the above approaches addresses the issue of quantum grav- 
itational loops. Since the loops can only be calculated within a well-defined per- 
turbation theory, one needs to construct a new perturbative expansion that would 
make diagrams tractable at short distances. 

Below we will rearrange perturbation theory in such a way that the consistent 
answers be obtained in the weak-field approximation. This can be achieved if the 
linearized gauge-fixing terms can play the role similar to the nonlinear terms. We 
will see that this requires a certain nontrivial modification of the linearized theory 
and of gauge-fixing procedure. 

We recall that in the DGP model the boundary (the brane) breaks explicitly 
translational invariance in the y direction, as well as the rotational symmetry that 
involves the y coordinate. However, this fact is not reflected in the linearized approxi- 
mation - the linearized theory that follow from (jHl) is invariant under five- dimensional 
reparametrizations.^^ This line of arguments suggests to introduce constraints in the 
linearized theory that would account for the broken symmetries. It is clear that an 
arbitrary set of such constraint cannot be consistent with equations of motion with 
boundary conditions on the brane and at — > oo. However, by trial and error a 
consistent set of constraints and gauge conditions can be found. Below we introduce 
this set of equations step by step. We start by imposing the following condition: 

= d^h, + d'^K^ = 0. (33) 

Furthermore, to make the kinetic term for the {/i5} component invertible we set a 
second condition: 

B, = d^h,, = Q. (34) 

At a first sight, the two conditions (jHHjl and (jH^ fix all the x-dependent gauge trans- 
formations and make the gauge kinetic terms nonsingular and invertible. However, 
at a closer inspection this does not appear to be satisfactory. One can look at the 
{yuz/} component of the equations of motion and integrate this equation with respect 
to y from — e to e, with e ^ 0. After the integration, all the terms with 5^ and B^ 
vanish. The resulting equation (which is just the Israel junction condition) taken 
by its own, is invariant under the following four-dimensional transformations 

h'^u{^\y)\y=o = hf,^{x,y)\y=o + df,Cu\y=0 + d^Cf,\y=o- (35) 

This suggests that in the rric —>■ limit the gauge kinetic term on the brane is 
not invertible. As a result, the problem of a precocious breakdown of perturbation 
theory discussed in the previous section arises. To avoid this difficulty one can 
introduce the following term on the brane world-volume: 

AS = -Ml,j d^xdy6{y) i^d^h^, - ^d^K^' . (36) 

^^If instead of the boundary we consider a dynamical brane of a nonzero tension, then the five- 
dimensional Poincare symmetry is nonlinearly realized and one has to include a Nambu-Goldstone 
mode on the brane. 
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This makes the graviton kinetic term of the brane invertible even in the 
hmit. At this stage, the partition function can be defined as 



Zgf = hma^^^o J dhAB exp (iS + i AS 



+ ,M?/dSd,{|i + |)). (37) 

Here S and AS are given in (jHJ and (jHUjl respectively, and the hmit a, 7 — > enforces 
flHHj) and (jnil). Before proceeding further, notice that Eqs. (jS^ and (jH^ would have 
been just gauge-fixing conditions if the boundary were absent (e.g., in a pure 5D 
theory with no brane). However, in the present case, the above equations, when 
combined with the junction condition across the brane, enforce certain boundary 
conditions on the brane. Therefore, Eqs. (jHHj) and ()34|) do more than gauge-fixing, 
and 7 and a cannot be regarded as gauge fixing parameters. The prescription given 
by ()37p is to calculate first all Green's functions and then take the limit a,7 0. 
Because of this, the results of the present calculations differ from ^3] where other 
boundary conditions were implied. 

Using (jH7|) we calculate below the propagator D and the amplitude A defined in 
fllO|) . We will see that there are no terms in D that blow up as rric —>■ 0. 

We start with the equations of motion that follow from ^7\. The {^i^} equation 
on the brane reads 



2 



r+e . 

J dy {dl,h^y - r]^^dl,h2 + d^d^h^ + d^d^h^^ - 2 r^^^S^^s/is 
T^u , (3J 



where 



In (I38|l we retained only terms that do not vanish in the e — >■ limit. Furthermore, 
G'^^}, denotes the 4D Einstein tensor, 

- Vf.udlhZ + v.ud^dph'^r (39) 
The {/iz^} equation in the bulk takes the form 

. (40) 
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At the next step we turn to the {/i5} equation which can be written as 

dlh^5 - d.d'^K^ - d.id'^h^, - d,hl) - - {d.d'^K^) = . (41) 

7 

Finally, the {55} equation takes the form 

dl K - d.d.h^^" - -{dl + d^d.hn = . (42) 
a 

The limit q;,7 — *• should be taken after the calculation is carried out. 

We turn to the momentum space with respect to four world- volume coordinates, 

hAB{x, y) = j d^pe'P^ JiABip, y) ■ (43) 

From the above equations we calculate the response of gravity to the source T^^. In 
the limit a, 7 — >■ the results are 

y) - -r^ — -\v,.T 4^1^) e-^\y\ . (44) 

We note that this expression is regular in the rric limit. This is contrary to 
what happens in the harmonic gauge Pl] where singular terms are present. 
For the off-diagonal components we find that ha5 ~ IPa, and 

ha5{p,y)^0. (45) 

Finally, 

- - ^/^^ - ^ ^-^ e-^l^l , (46) 

2 2 p^ + SrrLcP 



with r = (p^ + 2mcp)/(p^ + nicp)- The amplitude on the brane takes the form 

^l-graviton(P,2/ = 0) = ^ T - - T „ , ^ . (47) 



A remarkable property of this amplitude is that it interpolates between the 4D 
behavior at p ^ 



AMp.y = o) ^ - t;, - -t^) , (4j 



p2 V 2 

and the 5D amplitude at p -C mc 
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Branch cut 



Figure 2: The complex plane of Minkowskian momentum square with the branch cut and 
poles on a second Riemann sheet. The physical Riemann sheet is pole free. 



This amplitude has no van Dam-Veltman-Zakharov (vDVZ) discontinuity |^ 1^ 



It is instructive to rewrite the amplitude ()47|) in the following form: 



A- 



n 



graviton 



/2 



1 



+ rricP 6 p^ + TTicP 



where 



and 



T2 
1/2 



t;^ T ■ T 



SrUcP 



p"^ + SnicP 



(50) 



(51) 



(52) 



The first term on the right-hand side of (jSUj) is due to two transverse polarizations 
of the graviton, while the second term is due to an extra scalar polarization. The 
scalar acquires a momentum-dependent form-factor. The form-factor is such that 
at subhorizon distances, i.e., when p ^ rric, the scalar decouples. At these scales 
the effects of the extra polarization is suppressed by a factor rric/p (e.g., in the 
Solar system this is less than 10^^^). However, the scalar polarization kicks in at 
superhorizon scales, p <^ rric, where the five dimensional laws or gravity are restored. 



Let us discuss the above results in more detail. For this we study the pole 
structure of the amplitude (|5(J|) . There are two nontrivial poles 



p 



-rricP, 



and 



p 



—3 rricP . 



(53) 



19 



Let us find the positions of these poles on a complex plane of the Minkowski mo- 
mentum square p^, where = p^exp(— -ivr). For this we note that there is a branch 
cut from zero to plus infinity on the complex plane (see Fig. 2). The pole at = 
is just the origin of the branch cut with zero residue. Because of the cut the complex 
plane has many sheets (the propagator is multivalued function due to the square 
root in it). It is straightforward to show that both of the poles in fl53|) are on the 
unphysical, second Riemann sheet. Moreover, the positions of these poles are far 
away from the branch cut (usual particle physics resonances appear on unphysical 
sheets close to the branch cut, the above poles, however, are located on a negative 
semiaxis of the second Riemann sheet). Hence, the physical Riemann sheet is pole 
free. The poles on a nonphysical sheet correspond to metastable states that do not 
appear as in and out states in the S-matrix Using the contour of Fig. 2 that 
encloses the plane with no poles, and taking into account the jump across the cut, 
the four-dimensional Kallen-Lehmann representation can be written for the ampli- 
tude The latter warrants four-dimensional analyticity, causality and unitarity 
of the amplitude (j37j). Although the above interpretation is the only correct one, 
one could certainly adopt the following provisional picture that might be convenient 
for intuitive thinking. The second pole in (j53p can be interpreted as a "metastable 
ghost" with a momentum-dependent decay width that accompanies the fifth polar- 
ization and cancels its contributions at short distances. Remarkably, this state does 
not give rise to the usual instabilities because it can only appear in intermediate 
states in Feynman diagrams, but does not appear in the in and out states in the 
S-matrix elements. In this respect, it is more appropriate to think that the scalar 
graviton polarization acquires the form-factor g{p) (jK^ . 

The above results seem somewhat puzzling from the point of view of the Kaluza- 
Klein (KK) decomposition. Conventional intuition would suggest that the spectrum 
of the KK modes consists of massive spin-2 states. The Kallen-Lehmann represen- 
tation for the amplitude as a sum with respect to these massive states would give 
rise to the tensorial structure where the first term on the right-hand side of (j5U|) is 
proportional to Ty^, instead of Ty2- In this case, the remaining part of the am- 
plitude on the right-hand side would have a negative sign. This might be thought 
of as a problem. However, this is not so. The crucial difference of the present ap- 
proach from the conventional KK theories is that the effective 4D states are mixed 
states of an infinite number of tensor and scalar modes. What is responsible for the 
mixing between the different spin states is the brane-induced term and the present 
procedure of imposing the constraints. In the covariant gauge that we discuss the 
trace of h propagates and mixes with tensor fields. From the KK point of view this 
would look as an infinite tower of states with wrong kinetic terms. However, at 
least in the linearized approximation, the trace is a gauge artifact. Nevertheless, the 
effect of the trace part is that the true physical eigenmodes do not carry a definite 
four-dimensional spin of a local four-dimensional theory (see also |13]). Because of 
this there is no reason to split the amplitude ()50|) into the term that is proportional 
to T?/n and the rest. 
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The question of interactions of these states in the full nonlinear theory is not 
yet understood. What happens with the diagrams in which the "metastable ghosts" 
propagate in the loops (the unitarity cuts of which should give production of these 
multiple states) remains unknown. However, since the theory possesses 4D reparametriza- 
tion invariance, we expect that these questions will find answers similar to those of 
non-Abelian gauge fields. Further studies are being conducted to understand these 
issues. 

To summarize briefly, in this approach perturbation theory is well-formulated. 
The resulting amplitude interpolates between the 4D behavior at observable dis- 
tances and 5D behavior at superhorizon scales. This is due to the scalar polarization 
of the graviton that acquires a momentum-dependent form-factor. As a result, the 
scalar decouples with high accuracy from the observables at subhorizon distances. 

The model can potentially evade the no-go theorem for massive/metastable grav- 
ity that states that for the cancellation of the extra scalar polarization one 
should introduce a ghost that would give rise to instabilities jlHl d] • In the present 
model, at least in the linearized approximation, such instabilities do not occur. The 
convenient (although not precise) picture is to think of a "metastable ghost" that 
exists only as an intermediate state in Feynman diagrams which does not appear in 
the final states at least in the linearized theory. Since this state cannot be emitted 
in physical processes, it does not give rise to the usual instability. The latter prop- 
erty is similar to the observation made in the "dielectric" regularization of the DGP 
model in 

The questions that remain open concern the gauge-fixing and interactions in the 
full nonlinear theory where the Faddeev-Popov ghosts are expected to play a crucial 
role. These issues will be addressed elsewhere. 

3.3 Magic of nonlinear dynamics 

Exact static solutions in models of gravity carry a great deal of information on 
the gravitational theories themselves. Hence, finding these solutions in models that 
modify gravity at large distances is an important and interesting task. In this 
section, following Ref . ^] , we will study the Schwarzschild solution in the 5D DGP 
model jEj. It is comphcated to find this solution since even at distances much larger 
than the Schwarzschild radius of the source, full nonlinear treatment is required [lU] . 
The first approximate solution was obtained in Ref. and subsequently by the 
authors of Refs. jlHl EH EHl EH]- The solution should interpolate between very 
different distance scales. These scales are: the 4D gravitational radius of the source 
of a mass M, 

tm = 2GnM , 

the large distance crossover scale Vc ~ 10^*^ cm, and an intermediate scale, first 
discovered by Vainshtein in massive gravity [32], which in the DGP model reduces 
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HOI to 



n = [rMrl 



) 



1/3 



(54) 



This is a scale at which nonhnear interactions in a naive perturbative expansion 
in Gtv become comparable with the linear terms. For a source such as the Sun, 
the hierarchy of the scales is given in (j3U|) . Below, unless stated otherwise, we will 
consider sources smaller than r^. In Refs. [1711131131201123 approximate solutions 
for such sources were found in different regions of (jHUj) . The main properties of the 
solution can be summarized as follows: 

(a) At distances r ^ the 5D Schwarzschild solution with the 5D ADM mass 
M is recovered (throughout this work r stands for a 4D radius). 

(b) For r^. <C r <C Tc the potential scales as in the 4D Schwarzschild solution. 
However, relativistic gravity is a tensor-scalar theory that contains the gravitation- 
ally coupled scalar mode (i.e. the tensorial structure is that of a 5D gravitational 
theory which contains extra polarizations). 

(c) For r -C r* the theory reproduces the Schwarzschild solution of 4D general 
relativity (GR) with a good accuracy. 

Perhaps the most important property of the (a-c) solution outlined above is 
the dynamical "selfshielding" mechanism by which the solution protects itself from 
the would-be strong coupling regime pTOj. Very briefly, the selfshielding can be 
described as follows: the expansion in Gjv breaks down at the scale r ~ making 
the perturbative calculations unreliable below this scale. However, exact nonlinear 
solutions of equations of motion - which effectively resum the series of classical 
nonlinear graphs - are perfectly sensible well below the scale r^,. Hence, the correct 
way of doing the perturbative calculations is first to find a classical background 
solution of equations of motion and then expand around it. 

In Ref. jin] a 4D part of the metric was exactly found. This exact result, com- 
bined with reasonable boundary conditions in the bulk, is sufficient to determine 
unambiguously a number of crucial properties of the solution. First, this result 
confirms the existence of the scale r^^ - this scale enters manifestly our exact solu- 
tion. It also confirmed that the selfshielding mechanism outlined above takes place. 
Furthermore, it was emphasized that the selfshielding effect takes place because a 
source creates a nonzero scalar curvature that extends outside the source to a dis- 
tance r ~ r^.. This curvature suppresses nonlinear interactions that otherwise would 
become strong at the scale below r*. On the other hand, we also find that some 
of the physical properties of our solution differ from those in (a-c). The solution of 
Ref. [46] has the following main features: 

(A) For r ^ r^, like in (a), one recovers the 5D Schwarzschild solution, however 
unlike in (a), the new solution has the screened 5D ADM mass 




1/3 



(55) 
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The screened mass is suppressed compared to the bare mass M. Therefore, the new 
solution is energetically favorable over the (a-c) solution. 

(B) For -C r -C Tc one can think of the solution as being a four- dimensional 
one with an r-dependent decreasing mass M(r) ~ r^ruj'r- Alternatively, one can 
simply think of the solution just approaching very fast the 5D Schwarzschild metric 
with the screened mass fl55p . i.e., approaching the asymptotic behavior of (A). 

(C) For r <C r^, the results agree with those of (c) with a good accuracy. 

The (a-c) and (A-C) solutions both asymptote to the Minkowski space at infinity. 
However, the way they approach the flat space is different because of the difference in 
their 5D ADM masses. The (A-C) solution, or any of its parts, cannot be obtained 
in the linearized theory - it is a nonperturbative solution at any distance scale. 
Since the mass of the (a-c) solution is larger than the mass of the (A-C) solution, 
we would expect that the heavier solution will eventually decay into the light one, 
unless topological arguments prevent this decay. 

The above findings suggest that the Minkowski space, although globally stable 
in the DGP model, is locally unstable in the following sense. A static source placed 
on an empty brane creates a nonzero scalar curvature around it. For a source of the 
size ^ this curvature extends to a distance ~ r^. Above this scale the solution 
asymptotes very quickly to a 5D Minkowski space. More intuitively, a static source 
distorts a brane medium around it creating a potential well, and the distortion 
extends to a distance r ~ r*. Since is much bigger than the size of the source 
itself, we can interpret this phenomenon as a local instability of the fiat space. This 
local instability, however, has not been seen in the linearized theory [T^. It should 
emerge, therefore, in nonlinear interactions and should disappear when the scale 
tends to zero.^^ 

It is remarkable that the distance scale to which the local instability extends, 
coincides with the scale at which the naive perturbative expansion in Gn breaks 
down. Therefore, by creating a scalar curvature that extends to r ~ r*, the source 
shields itself from a would-be strong coupling regime that could otherwise appear 
at distances r <^ r^, [40 : (i) The coupling of a phenomenologically dangerous extra 
scalar polarization of a 5D graviton to 4D matter gets suppressed at distances r ^ 

due to the curvature effects. This is similar to the suppression of the extra 
polarization of a massive graviton in the AdS background jHIl E2] • Indeed, in our 
case the curvature created by the source, although coordinate dependent, has the 
definite sign that coincides with the sign of the AdS curvature. As a result, the model 
approximates with a high accuracy the Einstein gravity at r <^ r* with potentially 
observable small deviations EH] (see comments below), (ii) The selfcoupling 
of the extra polarizations of a graviton, which on a flat background leads to the 
breakdown of a perturbative expansion and to the strong couphng problem, gets 
now suppressed at distances ^ by the scalar curvature created by the source. 
This is also similar to the suppression of the selfcoupling of the massive graviton 

^■^The latter assertion is valid since the (A-C) solution, as we will see, is regular in the rric — > 
limit where it turns into a conventional 4D solution, i.e. it has no vDVZ discontinuity |48l I44| . 
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polarizations on the AdS background Elj • 

Above, we were primarily concerned with classical sources. Nevertheless, we 
would like to comment as well on dynamics of "quantum" sources, such as gravitons. 
Consider the following academic setup: a toy world in which there is no matter, 
radiation and/or any classical sources of gravity - only gravitons propagate and 
interact with each other in this world. Because of the very same trilinear vertex 
diagram that leads to the breakdown of the expansion for classical sources 
(see Ref. [101 ), the selfinteractions of gravitons will become important at lower 
energy scale than they would in the Einstein theory. The corresponding breakdown 
scale is the scale adopted to a quantum source with r^/ = 1/Mpi, that is 
A~^ ~ (rl/MpiY^^ [33j (see also Ref. that obtains a somewhat different scale). 
In this setup the graviton loop diagrams could in principle generate higher-derivative 
operators that are suppressed by the low scale. A theory with such high-derivative 
operators would not be predictive at distances below A^^ ~ 1000 km or so. 

However, there are two sets of arguments suggesting that the above difficulty 
might well be unimportant for the description of a real world which, on top of the 
gravitons, is inhabited by planets, stars, galaxies etc. We start with the arguments of 
Ref. [Sni- This work takes a point of view that A, is a true ultraviolet (UV) cutoff 
of the theory in a sense that at this scale some new quantum gravity degrees of 
freedom should be introduced in the model. Nevertheless, as was discussed in detail 
in Ref. |^, this should not be dangerous if one considers a realistic setup in which 
mater is introduced into the theory. For instance, consider the effect of introducing 
the classical gravitational field of the Sun. Because of the gravitational background 
of the Sun, the UV cutoff of the theory becomes a coordinate dependent quantity 
Ag{x). This cutoff grows closer to the source where its gravitational field becomes 
more and more pronounced, hence, increasing the value of the effective UV cutoff. 
In this approach the authors of Ref. managed to find a minimal required set 
of higher-dimensional operators that are closed with respect to the renormalization 
group fiow. Because of the resummation of large classical nonlinear effects these 
operators are effectively suppressed by the coordinate-dependent scale Ag{x). If so, 
the new UV physics will not manifest itself in any measurements [53] . 

Putting all this on a bit more general ground, one should define the model in 
an external background field. That is, in the action and the partition function of 
the model the metric splits into two parts (?^,^ = g'^^^^ + g"^^, where g"^^^ stands for 
the classical background metric and g'^^ denotes the quantum fiuctuations about 
that metric. The classical part satisfies the classical equations of motion with given 
classical gravitational sources such as planets, stars, galaxies etc.. Then, the effective 
UV cutoff for quantum fluctuations at any given point in space-time is a function 
of the background metric. For a realistic setup this effective cutoff is high enough 
to render the model consistent with observations. 

We flnd the above logic useful and viable. We also think that the algorithm of 
Ref. jSSj might be the most convenient one for practical calculations. Nevertheless, 
there could exist deeper dynamical phenomena beyond the above approach to the 
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discussions of which we turn right now. Although our arguments below parallel in a 
certain respect those of Ref. there is a conceptual difference on the main issue. 
Our view, that we will try to substantiate in subsequent works, is that the scale Ag 
is not a UV scale of the model in the sense that some new quantum gravity degrees 
of freedom should be entering at that scale. We think that all what's needed to go 
above the scale Ag is already in the model, and that this is just a matter of technical 
difficulty of nonperturbative calculations (or, in other words, is a matter of difficulty 
of summing up loop diagrams). The resummation could in principle cure problems 
at the loop level as well. At this end, we do not see a reason why the selfshielding 
mechanism outlined above should not be operative for "quantum" sources too. The 
very same local instability of the Minkowski space should manifest itself in nonlinear 
interactions of quantum sources, e.g., gravitons. The local instability scale in that 
case is Ag. Hence, we would expect that a quantum source creates a curvature 
around it that extends to the distances of the order of ~ (r^ 

/Mpi)i/3 ^ 1000 km, 

and doing so it selfshields itself from the strong coupling regime. If this is so, then the 
problem of loop calculations boils down to the problem of defining correct variables 
with respect to which the perturbative expansion should be performed. In this case 
the field decomposition should take the form: g^^, = g^]^ + ^f^^, where g'j^^ stands for 
a nonperturbative background metric created by a "quantum" source. Similar in 
spirit arguments using a toy model were given by Dvali in Ref. |56j . 

We find it useful to adopt a gauge in which the line element has an off-diagonal 
form: 

ds^ = e-^dt^ - e^dr^ - r^dQ^ - 7 drdy - e^dy^ , (56) 

where A, 7, a are functions of r and y. Our brane is located at ?/ = in this 
coordinate system. The Z2 symmetry w.r.t the brane implies that 7 is an odd 
function of y while the rest are even. A more conventional diagonal coordinate 
system can be obtained by a coordinate redefinition after which the interval reads 

ds^ = e'dt^ - e^dr^ - r^dn^ - e^dz" . (57) 

Here the functions u and P are related to A, 7, and a. In the z, r coordinate system 
the brane is bent. Typically in the brane- world models the 4D part of the Einstein 
equations are not closed. Hence, the induced metric on a brane cannot be deter- 
mined without some input from the bulk equations, and/or without making certain 
assumptions about the induced metric itself. This would also be true in our case. 
However, in the gauge (j56p . we find a subset of the Einstein equations that can be 
closed for the function A. As a result, A can be found exactly on the brane. Although 
the knowledge of A alone is not enough to describe all gravitational dynamics on the 
brane ( for instance, this is not enough for the description of the matter geodesies 
at short distances since transverse derivatives of the metric are also entering the 5D 
geodesic equations) nevertheless, combining the knowledge of A with the asymptotic 
behavior of the other functions in (|56p that we can also obtain unambiguously, is 
enough to deduce the properties (A-C). Hence, these properties are "exact." 
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Finally, we would like to make two important comments. First, the DGP model 
possesses two branches of solutions that are distinguished from each other by the 
bulk boundary conditions. These two branches are disconnected. In this work we 
concentrate primarily on the Schwarzschild solution of the so called conventional 
branch on which the brane and the bulk asymptote to the Minkowski space at 
infinity. However, the second, the so called "selfaccelerated" branch ^29j is extremely 
interesting as it can be used to describe the accelerated expansion of the Universe 
without introducing dark energy jHUj. In the present work we also find an exact 
brane metric for a Schwarzschild source on the selfaccelerated branch. However, 
because the asymptotic behavior of the solution on this branch is not Minkowski 
we are not able to argue for the existence of a nonsingular bulk solution. On the 
other hand, we do not see any physical reason why this solution should not exist 
in the bulk as well. This branch will be discussed in detail elsewhere. Second, it 
is interesting to note that the linearized analysis of the DGP model in dimensions 
six and higher jSTj, as well as certain modifications of the five- dimensional model 
[131 E] show no sign of breakdown of perturbation theory and strong nonlinear 
effects. It is left for future work to understand more deeply the interconnections 
between all these approaches. 

3.3.1 Structure of the solution 

In this subsection we discuss the properties of the solutions on the brane, i.e., at 
y = 0. In this discussion we closely follow Ref. [IH]. We find certain similarities, 
as well as drastic differences, in the 4D part of our solution with the anti-de Sitter- 
Schwarzschild (AdSS) solution of conventional 4D General Relativity (GR) with a 
small positive cosmological constant A (this is in spite of the fact that a source 
in our solution creates a curvature that has a signature of a negative cosmological 
constant) 



It is instructive to contrast our solution to the AdSS metric. 

Consider 4D GR with the cosmological constant A = — 3m^. Furthermore, con- 
sider a static source of mass M (a star) and a Schwarzschild radius tm = ^G^M in 
this space. In the static coordinate system the AdSS solution takes the form 



8 TT Gn T^y . 



(58) 





(59) 




This coordinate system covers the AdSS solution in the interval 



tm < f < Tc = 



-1 
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The following properties of the AdSS solution will be contrasted to our solution. 

(i) In the interval tm < r < there is a new distance scale r^, (|29|) exhibited 
by (|59|). The physical meaning of this scale is as follows. For r < r^, the Newtonian 
potential tm/t in ()59|) dominates over the term mlr"^, while for r > r^, the term 
m^r^ overcomes the Newtonian term. Hence, r^, is a scale at which the Newtonian 
and the m^r^ terms are equal. This can also be expressed in terms of invariants. 
Let us define the Kretschman scalar (KS) 



Rk ^ ^J{R^%sY , (60) 

where R^'fj^yg is a Riemann tensor of the Schwarzschild part of the solution (i.e., of the 
part that survives in the limit). We compare the KS with the background 

curvature due to the cosmological constant 

\Ra\ = 12ml. (61) 

We get 

Rk > \Ra\ for r < r, ; i?^ < \Ra\ for < r < . (62) 

Therefore, r* is a scale at which Rk — \Ra\- For tm r r^, the corrections 
due to the background curvature are small and the solution is dominated by the 
Schwarzschild metric, while for r^, ^ r ^ the background curvature terms are 
larger that the Schwarzschild terms, both of them still being smaller than 1. 

(ii) At r > Tc the Schwarzschild part becomes irrelevant compared to the AdSS 
part. 

We will show below that our solution has some of the properties described in (i), 
however, unlike (ii), it behaves as 5D Schwarzschild solution at large distances. 
The 4D part of our solution (i.e. the solution at ?/ = 0) for r <^ takes the form 

'^''^ -r^dnj. (63) 



(1 



rM 



Like the AdSS solution, the metric (jMjl possesses the scale defined in (jSHI)- As 
we will see below this scale has the same physical meaning as in the AdSS case. For 
instance, at r <^ r* 

9{r) ^ [f,^ . (64) 

Then, it is straightforward to check that 

— ^ ^ mlr'^gir) for r <^ r*; — ^ ~ mlr'^g{r) for r ~ . (65) 
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Figure 3: The solid line depicts P{r)/rM (on the vertical axis, where P = r(l — (70O)), 
as a function of r on the horizontal axis. The dashed line presents the function 0.28r*/r; 
The value of r* is set to 1 on this graph. 

Therefore, the corrections become of the order of the Tm/t term at around r ~ r^,. 
Moreover, like the AdSS solution, the corrections dominate over Tm/t for r^, ^ r ^ 
Tc turning the 4D behavior of the solution into the 5D behavior. The plot of the 
function is given on Fig. 3. As in the AdSS case, the corrections to the Schwarzschild 
solution that are proportional to rUc give rise to the four-dimensional Ricci curvature 
Rm^. This is interesting since the curvature is completely due to the modification 
of gravity. However, unlike the AdSS case, this curvature is not a constant but 
depends on r; moreover it also depends on the strength of the source itself. The plot 
of the Ricci curvature is given on Fig. 4. The presence of this curvature can easily 
be understood by looking at the trace of the 4D Einstein equation on the brane 

i?4 + ?>mcK = SttGnT. (66) 

T is zero outside a localized source such as a star. However, the trace of the extrinsic 
curvature K is not zero, therefore, = — 3mc K ^ and i?4 outside of the source 
is nonzero and equals to Rm^. 

Similar to the AdSS solution the above properties can be expressed in terms of 
the invariants 

Rk > Rm^ for r < and r > ; Rk ^ Rmc for r . (67) 

Unlike the AdSS solution, however, the curvature Rm^(r) decreases very fast after 
r 3> r*. Hence, the induced curvature Rm^{r) is subdominant to Rk everywhere 
except in the neighborhood of the point r ~ where they both are of the same 
order ~ m^, see Fig. 4. 

Furthermore, unlike the AdSS solution, our solution can be presented in the same 
coordinate system even for r > Tc- This is because there is no horizon at r = Tc and 
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Figure 4: The solid line depicts the magnitude of the four-dimensional Ricci scalar cur- 
vature (on the vertical axes) as a function of r on the horizontal axes. The dashed line 
depicts the dependence of the 4D Kretschmann scalar on r. The value of r* is set to 1 on 
this plot. 



our solution smoothly turns into the 5D Schwarzschild solution 

ds\=o = fl - %1 - - r^d^l. (68) 

The key property of this solution is that the gravitational radius is rescaled 

tm ^ tm [ — > r-M • (69) 

This has an explanation. The gravitational radius grows compared to tm because 
in the 5D regime the gravitational coupling constant grows. However, there is an 
opposite effect as well. In fact, the gravitational radius reduces compared to what 
it would have been in a pure 5D theory with no brane. This is because the effective 
mass of the source M^g, defined as rj^ = Mes/M^, gets screened. Indeed, M^s 
includes contributions from the curvature Rm^{r) that stretches out all the way to 
r — rc- The screened mass of the source is 

Meff ~ M f!JiY^ < M. (70) 

Thus, as seen from r ^ Tc distances, there is strong screening of the source. 

All the above results could be understood as follows. Consider an empty brane 
and an empty bulk. The Minkowski space is a solution. Let us place a static source 
on the brane at r = 0. The Minkowski solution remains globally stable, however, 
the source, no matter how weak, triggers local instability of the Minkowski space on 



29 



a brane in the region r < r*. In this patch the Minkowski space is readjusted to a 
curved space. The curvature of the latter depends on the strength of the source, it 
slowly decreases with increasing r but drops fast at r > r*. For an observer at large 
distances it looks as if the source has polarized the medium (brane) around it. This 
observer measures the screened mass (jU^ which also includes the contributions of 
the curvature. 

At large enough distances, i.e. i/r^ + ^ r^, the solution turns into a 5D 
spherically-symmetric Schwarzschild solution. 



However, the 5D spherical symmetric is only an approximation and does not hold 
for \/r^~+lj^ <^ r^. In the latter regime the properties of the solution on and off 
the brane are rather different. The pure 5D spherically-symmetric solution (f7T|l is 
squeezed both on and off the brane but it is more squeezed on the brane than in 
the bulk. Hence, the only symmetry of the solution that is left is the cylindrical 
symmetry. 

4 Brane-induced gravity in more than five dimen- 
sions — softly massive gravity 

In the present section we turn to Brane Induced Gravity in more than 5 dimensions. 
These models share the main property of the 5D model that gravity at short dis- 
tances is four-dimensional, becoming higher dimensional at larger scales. However, 
there are a number of crucial distinctions from the 5D case. These are: 

• Unlike a 3-brane embedded in 5D space-time, the embedding into six and 
higher dimensional space-times allows to preserve the 3-brane world-volume 
to be flat. This suggest that the CCP could be solved in this approach in 
dimensions six and higher fT^ 1^ . 

• Unlike the 5D case, the transverse to the brane Green's functions in six di- 
mensions and higher are sensitive to the UV physics. A careful treatment of 
the UV regularization procedure [13 EHl IMl 1^0] is needed. The best option is 
to use the UV completion dictated by the String Theory construction 26 . 

• In most of the interesting cases where simplest UV regularizations were tried so 
far, the fiat space propagator exhibits new poles that correspond to very light 
m ~ i^o tachyons WD^ with negative residues [HI] ■ The positions of these poles 






(71) 
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and their mere existence is UV regularization dependent. At present, it is not 
clear whether these poles would persist in a consistent UV completed theory, 
such as in the string theory construction j^H]- Work to establish this is in 
progress [021 • However, even in the regularizations where the new poles exist, 
they can be eliminated by a special prescription for the poles of the Greens 
functions ^\ at the expense of sacrificing causality at the Hubble scales. This 
suggest that if these poles are truly present in the theory, then the flat space- 
time, although globally stable, might not be a locally stable background within 
a patch of the size Hq^. Such a local instability takes place in the 5D model 
but only at a nonlinear level and at a different scale jlHl- This issue in the 
context of six- and higher- dimensional models awaits further investigation. 

• Unlike the 5D case, it turns out that in six and higher dimensions perturbation 
theory does not break down at low scale [57^ . This observation is independent 
of the UV regularization as well as of the presence/absence of the new poles 
discussed in the previous item (i.e., the naive perturbation theory does not 
break down at a low scale even when the new poles are not present, see below). 

The first issue on the list above was already discussed in Sect. 2. In what 
follows we will study in turn the second, third and fourth issues. As was mentioned 
above, these three items apply to a brane that has no tension (i.e, to a brane the 
cosmological constant on which is fine-tuned to zero). However, our main goal is to 
study a brane with an arbitrary tension, and the last three properties listed above 
should be studied and established for this case. This is only partially fulfilled so far, 
with certain encouraging results (see below). Further detailed calculations are still 
needed. Here we derive results for a tensionless brane following Ref. jS^. 

The equation of motion for the theory described by the action ^ in D > Q takes 
the form 

5(^)(l/)M^iG(f)<5!;<5^ + M^^gSS = -T,^{x)5>:,5^^S^''\y). (72) 

G^i^} and G^^^ denote the four-dimensional and Z)- dimensional Einstein tensors, re- 
spectively. We choose (for simplicity) a source locahzed on the brane, T^y{x)5'^^\y). 

Gravitational dynamics encoded in Eq. ()72|) can be inferred both from the four- 
dimensional (4D) as well as (4 + A^)-dimensional standpoints. From the 4D perspec- 
tive, gravity on the brane is mediated by an infinite number of the Kaluza-Klein 
(KK) modes that have no mass gap. Under conventional circumstances (i.e., with 
no brane kinetic term) this would lead to higher- dimensional interactions. However, 
the large 4D Einstein-Hilbert (EH) term suppresses the wave functions of heavier 
KK modes, so that in effect they do not participate in the gravitational interac- 
tions on the brane at observable distances JBT. Only light KK modes, with masses 
rriKK ^ "^c, 

Ml , , 
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remain essential, and they collectively act as an effective 4D graviton with a typical 
mass of the order of and a certain smaller width. 

Assuming that M^, ~ 10~^ eV or so, we obtain ~ Hq ~ 10^^^ GeV. Therefore, 
the DGP model with N >2 predicts [^01 a modification of gravity at short distances 
~ 0.1 mm and at large distances ~ Hq^ ~ 10^*^ cm, give or take an 
order of magnitude. Since gravitational interactions, nevertheless, are mediated by 
an infinite number of states at arbitrarily low energy scale, the effective theory (jH)) 
presents, from the 4D standpoint, a non/oca/ theory [201 ■ Moreover, as was suggested 
in [21], nonlocalities postulated in pure 4D theory can solve the CCP and give 
rise to new mechanisms for the present-day acceleration of the universe |2I]. (It is 
interesting to note that the nonlocalities in a gravitational theory that are needed 
to solve the cosmological constant problem could appear from quantum gravity jHSl 
or matter loops [Hni in a purely 4D context.) 

On the other hand, from the (4 + A^)-dimensional perspective, gravitational 
interactions are mediated by a single higher-dimensional graviton. This graviton 
has two kinetic terms given in Eq. (jH)), and, therefore, can propagate differently on 
and off the brane. Namely, at short distances, i.e. at r < ~ Hq^ ~ 10^^ cm, 
the graviton emitted along the brane essentially propagates along the brane and 
mediates 4D interactions. However, at larger distances, the extra-dimensional effects 
take over and gravity becomes (4 + A^)-dimensional. 

As was first argued in Ref. the results in iV > 2 DGP models are sensitive 
to ultraviolet (UV) physics, in contradistinction to the = 1 model (14j. In other 
words, one should either consistently smooth out the width of the brane jHOl, or 
introduce a manifest UV cutoff in the theory |HH1 tHHl lEHI , or do both. With a finite 
thickness, more localized operators appear on the world-volume of the brane, in 
addition to the world- volume Einstein-Hilbert term already present in Eq. (jH} jBUj . 
For instance, one could think of a higher- dimensional Ricci scalar smoothly spread 
over the world- volume j45j . 

In general, terms that are square of the extrinsic curvature can also emerge. 
Some of these terms can survive in the limit when the brane thickness tends to zero 
(i.e. in the low-energy approximation). For instance, in the zero-thickness limit of 
the brane the following terms might be important: 

6^''\y)h';dlh:, 6^''\y)h^''d,d,K, S^'^^y) h'^dad.h^' , (74) 

where h denotes small perturbations on flat space. Although the main features of the 
model, such as interpolation between the 4D power-law behavior of a nonrelativistic 
potential at short distances and the higher- dimensional behavior at large distances, 
are not expected to be changed by adding these terms, nevertheless, the tensorial 
structure of a propagator could in general depend on these terms and selfconsistency 
of the theory may require some of these terms to be present in the actions in a 
reparametrization invariant way. 

In the low-energy approximation the exact form of these "extra" terms and their 
coefficients are ambiguous, because of their UV origin. They will be fixed in a 
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fundamental theory from which the DGP model can be derived 1^ IHHj . In the 
present paper, in the absence of such a fundamental theory, (but in the anticipation 
of its advent), we would like to study a particular parametrization of these "extra" 
terms, for demonstrational purposes. According to our expectations, physics in 
the selfconsistent theory will have properties very similar to those discussed below. 
We will show that these properties are rather attractive since they do avoid severe 
problems of 4D massive gravity. 
Consider the action 

S = Ml, j d*x^{aR{g) + 67^4+^) 

+ M^+'' J d'xd^'y^n.Mg), (75) 

where, in addition to the 4D EH term, a D-dimensional EH term localized on the 
brane is present. Here a and b are some numerical coefficients. We will study the 
properties of the system described by ()75p for different values of a and b. The 
action (f7H|) is fully consistent with the philosophy of Ref. [13]: if there is a (1+3)- 
dimensional brane in D-dimensional space, with some "matter fields" confined to 
this brane, quantum loops of the confined matter will induce all possible structures 
consistent with the geometry of the problem, i.e. (l+3)-dimensional wall embedded 
in D-dimensional space. 

The equation of motion in the model ()75j) takes the form 

+ = -T,.(x) 6^^ 5), 5(^)(y) . (76) 

In deriving the above equation we first introduced a finite brane width A, and then 
took the A — »• limit in such a way that no surface terms appear. In general, the 
results depend on the regularization procedure for the brane width. In the present 
work we adopt a simple prescription in which derivatives with respect to the trans- 
verse coordinates calculated on the brane vanish in the A ^ limit (a unique 
prescription could only be specified by a fundamental theory.). As previously, G*-^-* 
and G'^^^ denote the four-dimensional and D-dimensional Einstein tensors, respec- 
tively, while a and b are certain constants. In order to be able to describe 4D gravity 
at short distances with the right value of Newton's coupling we set 

a + b = I. {11) 

Note, that the first two terms in parenthesis on the left-hand side of Eq. ()76p can 
be identically rewritten as 

(a + 6)Gg) +b(-d,daK - d^dahl - + v^.d^h^ 

+ d^d^K - r]^,^dl K + 2r]^,dadah'''' + r]^Mh''') . (78) 
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The above equation of motion (j76|) - which should be viewed as a regularized version 
of the DGP model - could be obtained from the action (jSI) as well, provided the 
latter is amended by certain extrinsic curvature terms. Below we will study this 
version of the regularized DGP model for various values of the parameters a and h. 



4.1 Perturbation theory of flat space 

We start our discussion with a simple model of a scalar field $ in (4 +A^) -dimensional 
space-time. For convenience we separate the dependence of the scalar field $ on 
four- dimensional and higher- dimensional coordinates, $(x^,ya) = ^{x^y)- The two 
kinetic term action - the scalar counterpart of Eq. (jS} - is 

S = M|i y"ci^a;a^$(a;,0)a'^$(x,0) 

+ Ml+^ J d^x d^y dA^{x, y) d^^{x, y) . (79) 

It is important to understand that in the scalar case the analog of the new term 
included in Eq. (fTUj) but absent in (f?^ reduces, identically, to the already existing 
localized term. This is a consequence of our choice of the regularization of the brane 
width A and the boundary conditions according to which transverse derivatives 
vanish on the brane in the A — > limit. 

To study interactions mediated by the scalar field we assume that $ couples to 
a source J localized in the 4D subspace in a conventional way, / (i^x2$(x, 0) J{x). 
Then the equation of motion takes the form 

(y) Mil a; $(x, 0) + M2+^ d\^{x,y) = J{x) {y) . (80) 

The very same equation applies to the scalar field Green's function. 

To solve this equation it is convenient to Fourier-transform it with respect to 
"our" four space-time coordinates x^ p^, keeping the extra y coordinates intact. 
Marking the Fourier-transformed quantities by the tilde, 

<^{x,y) <l(p,y), (81) 

we then get from Eq. (jSUI) 

5W(y)M|i(V)$(p,0) 

+ Mr^ (V - \mp.y) = -Hp) s^'^Ky) , (82) 

where = Po — pf — P2 " Ph ^''^^ the notation 

N Q2 

A» - E g-i (83) 

a=l ^i/a 
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is used. 

We will look for the solution of Eq. (j82j) in the following form: 

= D{p,y)x{p), (84) 
where the function D is defined as a solution of the equation 

^-p2_^^^_,^^Dip,y) = (85) 

Note that the function D is uniquely determined only after the ie prescription spec- 
ified above is implemented. We also introduce a convenient abbreviation 

Doip) = D{p,y = 0). (86) 



Now, it is quite obvious that a formal solution of Eq. ()82j) can be written in 
terms of the function D as follows: 

y) = - M ^(P^y^ + cl., (n y) (87) 

where ^horaij>^y) is ^ general solution of the corresponding homogeneous equation 
(i.e., Eq. (jH^ with the vanishing right-hand side), and c is an arbitrary constant. 
Equation (jHTj) presents, in fact, the Green's function too, up to the factor J(p)/Mpi, 
which must be amputated. In particular, for the Green's function on the brane we 
have 

G{p, 0) = ^ $(p, y = 0) , Ghom(p, 0) = ^ $hom(p, 1/ = 0) , (88) 

J [p] J [p) 

while for arbitrary values of y 



p"^ Do{p) — 



GiP,y) = - ^,r. )J^ \n + cGhom(p,y) , (89) 



where 



N 



u^^'-!^ = mlM^-\ (90) 

The presence/absence of the homogeneous part is regulated by the ie prescription. 
Note that if the first term on the right-hand side of Eq. (jH^ has poles on the real 
axis of p^, then the homogeneous equation has a solution 

G'hom(p, y) = D{p, y) 6 {fD,{p) - M^) . (91) 

This fact will play an essential role for gravity. 

In what follows we will examine the poles of the Green's function G{p,y). The 
positions of these poles depend on the functions (jhom(p, 2/), and Dq as defined in 
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Eqs. fl85p and ()86p. The choice of a particular rule of treatment of the poles cor- 
responds to the choice of appropriate boundary conditions in the coordinate space. 
Note that the latter are dictated by physical constraints on the Green's function G 
rather than on the auxiliary function D. 

To get to the main point, we will try the simplest strategy of specifying the poles 
and check, a posteriori, whether this strategy is self consistent. Let us put 

c = 

and define D in the Euclidean momentum space. Since in the Euclidean space the 
expression for D is well-defined and has no singularities, 

D{pE^q) = , D{pE,q) = / d'^y e'^y D{pE,y) , (92) 

Pe -r q J 

a 

one can perform analytic continuation from the Euclidean space to Minkowski. This 
is not the end of the story, however. It is the Green's function G that we are 
interested in, not the auxiliary function D. As will be explained below, the above 
procedure is consistent, for the following reason. The function G obtained in this 
way has a cut extending from zero to infinity. In addition, we find two complex 
conjugate poles on the second unphysical Riemann sheet of the complex p"^ plane. 
Moreover, there are additional poles on subsequent unphysical sheets. 



4.2 Six dimensions 

It is instructive to demonstrate how things work by considering separately the six- 
dimensional case. In six dimensions sensitivity to the UV cutoff is only logarithmic, 
and it is conceivable that the results obtained in the cut-off theory could be consis- 
tently matched to those of a more fundamental UV-completed theory-to-come. 
It is not difficult to calculate 

Do(s) =^ln(^^ +lj , s=p\ (93) 

where is an ultra-violet cut-off. With this expression for Dq the function G(p^, 0) 
develops a cut on the positive semiaxes of s due to the logarithmic behavior of Do{s). 
This fact has a physical interpretation. Since the extra dimensions are noncompact 
in the model under consideration, the spectrum of the theory, as seen from the 4D 
standpoint, consists of an infinite gapless tower of the KK modes. This generates a 
cut in the Green's function for s ranging from zero to +oo. 

^^The D > 6 models of brane-induced gravity are power sensitive to UV physics. In general one 
expects all sorts of higher-derivative operators in this case. 
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In addition, there might exist isolated singular points in G(p^,0). These singu- 
larities (for s <^ A^) are determined by the equation 



s — m,. 



A-n \-s , 



0. 



(94) 



where ml is defined in Eq. flTHj) . Let us introduce the notation 

s = Soexp(27) , (95) 
where Sq is a real positive number. Then, Eq. (jMj) has two solutions of the form 



So ~ 47rm 



1 

mi 



and 



TT 



7l 



log(A2/m^ 



72 ^ 27r 



vr 



log(A2/m 



(96) 



(97) 



We conclude that there are two complex-conjugate poles on the nearby unphysical 
Riemann sheets. These poles cannot be identified with any physical states of the 
theory. They are, in fact, manifestations of a massive resonance state. All other 
complex poles appear on subsequent nonphysical Riemann sheets. 



4.3 More than six dimensions 

Physics at D > 6 is similar to that of the six- dimensional world which was de- 
scribed in Sect. 14.21 There are minor technical differences between odd- and even- 
dimensional spaces, however, as we will discuss momentarily. 
In seven dimensions we find 



Do{s) = ^ |A - v^arctan (^y^^J | • (98) 

As in the 6D case, there is a branch cut. The cut in this case is due to the depen- 
dence of the Green's function on ^/s. No other singularities appear on the physical 
Riemann sheet. All poles are on unphysical Riemann sheets, as previously. 
In the eight-dimensional space the expression for Do reads 

Again, we find a cut due to the logarithm, similar to that of the 6D case. All isolated 
singularities appear on unphysical Riemann sheets. 
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The nine-dimensional formula runs parallel to that in seven dimensions, 



Do{s) = i ^ ~ v^arctan 





(100) 



Finally, in ten dimensions 



Do{s) = 1— + s\a^ + s fin— + . 

' 128 7r3 \ 2 [ V /J J 



(101) 



The pole structure of G is identical to that of the eight-dimensional case. Since the 
pattern is now well established and clear-cut, there seems to be no need in dwelling 
on higher dimensions. 

Before turning to gravitons we would like to make comments concerning the 
UV cutoff A. The crossover distance Vc ~ depends on this scale: in 6D the 
dependence is logarithmic, while in D > 6 this dependence presents a power law 
Hence, the crossover scale in the N >2 DGP models, unlike that in the 

= 1 model, is sensitive to particular details of the UV completion of the theory. 
Since in the present work we adopt an affective low-energy field-theory strategy, we 
are bound to follow the least favorable scenario in which the cutoff and the bulk 
gravity scale coincide with each other and both are equal to M^, ~ 10^^ eV. If a 
particular UV completion were available, it could well happen that the UV cutoff 
and bulk gravity scale were different from the above estimate. In fact, in the string- 
theory-based construction of Ref. the UV completion is such that the cutoff 
and bulk gravity scale are in the ballpark of TeV. 

In conclusion of this section it is worth noting that the Green's function Dq in the 

> 3 case contains terms responsible for branch cuts. These terms are suppressed 
by powers of s/A, and, naively, could have been neglected. It is true, though, that 
the explicit form of these terms is UV-sensitive and cannot be established without 
the knowledge of UV physics. One should be aware of these terms since they reflect 
underlying physics - the presence of the infinite tower of the KK states. Fortunately, 
none of the results of the present work depend on these terms. 

4.4 The graviton propagator 

Now it is time to turn to gravitons with their specific tensorial structure. We will 
consider and analyze the equation of motion of the DGP-type model presented in 
Eq. ()76|) . which we reproduce here again for convenience 



Here G'-^^ and G^^^ denote the four-dimensional and D-dimensional Einstein tensors, 
respectively, while a and b are certain constants satisfying the constraint 

a + b = l. 



+ Mr^Gi^i = -T,.(x)5^5^5(^)(y). 



(102) 
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For simplicity we choose a source term localized on the brane, namely, T^^{x)S^^\y). 
At the effective-theory level the ratio a/b = a/(l — a) is a free parameter. The 
only guidelines we have for its determination are (i) phenomenological viability; 
(ii) intrinsic selfconsistency of the effective theory which, by assumption, emerges 
as a low-energy limit of a self consistent UV-completed underlying "prototheory." 
Specifying the prototheory would allow one to fix the ratio a/(l — a) in terms of 
fundamental parameters. 

Our task is to study the gravitational field produced by the source T^^{x)5^^\y). 
To this end we linearize Eq. (|1U2|) . If g^B = Vab + '^hAs, in the linearized in h 
approximation we find 

gSJ = dlhAB - dAdch% - dndch'i 

+ Oa Ob h^c - VAB dl + tjab dc do h""^ , (103) 

where djy = dod^ . On the other hand, the four-dimensional Einstein tensor in the 
linearized approximation is 

- V,udlK + V^udc^dph^^. (104) 
In what follows we will work in the harmonic gauge, 

d^hAB = \dBh''c. (105) 

The advantage of this gauge is that in this gauge the expression for G^ab significantly 
simplifies, 

= dl hAB - \ VAB dl h'^l . (106) 

Additional conditions which are invoked to solve the {ab} and {a/i} components of 
the equations of motion are 

ha^, = 0, hab = ^Vabhc- (107) 
Using the last equation it is not difficult to obtain the relation 

Nh>;^ = {2-N)h:. (108) 

This relation obviously suggests that we should consider separately two cases: 

(i) N = 2; 

(ii) > 2. 

We will see, however, that the results in the N = 2 and N > 2 cases are somewhat 
similar. 
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4.5 Brane-induced gravity in six dimensions {N = 2) 

In two extra dimensions Eq. (jlOHj) implies 

h', = 0- (109) 

Therefore, the trace of the D- dimensional graviton coincides with the trace of the 
extra-dimensional part, 

hi = K. (110) 

As a result, the four-dimensional components of the harmonic gauge condition p05|) 
reduce to 

d^h,, = ^d,K. (Ill) 

Let us now have a closer look at the {fiiy} part of Eq. (fTBj) . Taking the trace of this 
equation and using Eqs. ()106|) . ()104|) . ()109|) and ()111|) we arrive at^^ 

(36 - l)5W(y)M|,9j/^: + 2M^^9>: = T;5(^)(y). (112) 

The obtained equation is very similar to the scalar- field equation (j80|) . Therefore, 
we will follow the same route as in the scalar-field case, until we come to a subtle 
point, a would-be obstacle, which was nonexistent in the scalar- field case. 
Let us Fourier-transform Eq. ()112|) . 

{3b-l)5^''\y)MU-p')K{p,y) 

+2M^^ (V - A,) Kip, y) = f{p) 5^(2/) . (113) 
The general solution of the above equation is 



K{p.y) = ^Q{p.y). im 



Tip} 

2 

D{p,y) 



+ C^hom, (115) 



2ml - (36 - l)p'^Do{p) 
where the solution of the homogeneous equation takes the form 

Ghom = Dip,y)5(2ml - (36 - l)//^o(p)) • (116) 



To begin with, let us consider the case 36 > 1. Then the first term on the right-hand 
side of Eq. flll5|l has poles for complex values of p"^. For instance, in the 6D case 
this pole is determined by the equation 



2mj An 2ml 



(36 - 1)Dq{s) (36 



—s 



;ii7) 



15 



As before, we put the transverse derivatives to be zero in the A ^ hmit. 
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This equation has at least two solutions of the form 



An 2ml 
3b- 1 



In- 



A2 



and 



TT 



71 



log(A^/m^ 



72 ^ 27T 



TT 



log(A^/m 



118) 



(119) 



The quantity h'^{p,y) is not a gauge invariant variable. Therefore, the presence of 
certain poles in the expression for h'^{p,y) depends on a gauge. However, exphcit 
calculations (see below) show that the poles found above also enter the gauge in- 
variant physical amplitude. Therefore, we need to take these poles seriously and 
analyze their physical consequences. 



4.6 6>l/3 

If 6 > 1/3 there are no poles on the physical Riemann sheet. Instead, poles appear 
on the nearest unphysical Riemann sheets. These poles cannot be identified with 
any physical states of the theory. They represent a signature of massive resonance 
states. All other complex poles appear on subsequent nonphysical Riemann sheets. 

Using a contour integral one can easily write down the spectral representation 
for the Green's function Q 

a(p,, = o) = ir-JiWi?L^, (120) 

TT JO t — p"^ — te 



where the spectral function is defined as 

[(36 - l)t Re Do - 2m^f + [(36 - l)t Im Dq 

and 



JV-2 



We see that p(t) satisfies the positivity requirement. Equation p2()j) guarantees that 
the Green's function Q is causal. 

The next step is applying the expression for Q to calculate h^j^. In fact, it is 
more convenient to calculate the tree-level amplitude 

A{p,y) ^ K,{p,y)f'>^%p), (123) 

where T'^'^[p) is a conserved energy momentum tensor, 

p^f'^"' = p^f'^'' = 0. 
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Using Eqs. ()102p . ()114|) and ()134|) we obtain a lengthy expression for the amphtude 



A(v v) - if f 



TV 



(26 - lyz^o - ml 
(36 - l)p2Do - 2ml 



(124) 



Let us study the above expression in some detaiL The first question to ask is about 
poles. It is quite clear that the p^-poles of A are of two types; their position is 
determined by: 



c 

or 



{3b - l)p'^ Do{p) = 2m; 



As was explained previously, all these poles appear on the second Riemann sheet, 
with the additional images on other unphysical sheets. None of these poles can 
be identified with asymptotic physical states. As was elucidated above, the occur- 
rence of the poles on the second and subsequent Riemann sheets corresponds to the 
massive-resonance nature of the effective 4D graviton. Our previous analysis can 
be repeated practically verbatim, with minor modifications, proving analyticity and 
causality of the amplitude A. 

Next, we observe that at large momenta, i.e., when p"^ Dq{p) ^ ml, the scalar 
part of the propagator has 4D behavior; the tensorial structure is not four-dimensional, 
however. The terms in the braces in Eq. ()124|) . namely. 



correspond to the exchange of massive gravitons and scalar degrees of freedom. This 
would give rise to additional contributions in the light bending, and is excluded 
phenomenologically, unless the contribution due to extra polarizations is canceled 
by some other interactions (such as, e.g., an additional repulsive vector exchange). 
Note also that when b ^ a, i.e., 6 — > 1, one obtains the tensorial structure of 6D 
gravity, as expected from (f7^ . 

On the other hand, at large distances, i.e., at -Do(p) ml,'we get the following 
tensorial structure of the amplitude (jl24j) : 

f^uf'^" - f'l . (126) 

This exactly corresponds to the exchange of a six-dimensional graviton, as was 
expected. 
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4.7 6<l/3 

This case is conceptually different from that of Sect. 14. 61 As we will see momentarily, 
if 6 < 1/3 there are no problems in (i) maintaining 4D unitarity; and (ii) getting 
the appropriate 4D tensorial structure of gravity at subhorizon distances. This is 
achieved at a price of abandoning 4D analyticity, in its standard form, which could 
presumably lead to the loss of causality at distances of the order of ~ 10^*^ cm. 
The absence of causality at distances > 10^^ cm, was argued recently to be an 
essential ingredient for solving the cosmological constant problem. 

Although all derivations and conclusions are quite similar for any ratio a/6 as 
long as 2b < a, we will stick to the technically simplest example 6 = 0, a = 1. In 
the situation at hand, the homogeneous part pifjj) need not be trivial, i.e c need 
not vanish. The value of the constant c is determined once the rules for the pole 
at p^Dq{p) + 2ml — specified. Putting c = leads to nonunitary Green's 
function. Therefore, we abandon the condition c = in an attempt to make a more 
consistent choice that would guarantee 4D unitarity. We stress that we are after 
unitarity here, not unitarity plus causality. 

To begin with we pass to the Euclidean space in (i.e. p\) and introduce 

the following notation: 

^^""^^^^^ " 2ml - plD.ip,,) - ^e ' ^^^^^ 

The function P^^'^ is a Euclidean-space solution of Eq. (|113|) . with the particular 
choice c = in. (The choice c = —in would lead to Eq. p27|) with the replacement 
e^-e). 

As the next step we will analyze the complex plane of p\. Since the function 
Dq{pe) is real, the function P^^\pe) must (and does) have an isolated singularity 
in the p'^ plane which is similar to a conventional massive pole, except that it lies 
in the Euchdean domain. This singularity occurs at the point p"^ = pi is defined by 
the condition 

pIDo{p^) = 2ml, pI real and positive. (128) 

This is the only isolated singularity in Eq. ()127p : it is located in the complex j9|; plane 
on the real positive semiaxis. In addition to this pole singularity, the function ()127|) 
has a branch cut stretching from zero to —oo due to the imaginary part of Dq{p\) 
appearing at negative values of p\. As before, this branch cut is the reflection 
of an infinite gap less tower of the KK states. As a result, the following spectral 
representation obviously emerges for P''^\p\): 



TT Jo U - p% Vi - Pe - 

with the Euclidean pole term being "unconventional." The residue of the pole R is 
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given (for any A^) by 



Note that in the first term on the right-hand side of Eq. p29p the integration 
runs from zero to minus infinity; thus, the integrand never hits the would-be pole at 
u = p\> Therefore, the ie prescription is in fact used only to specify the isolated 
pole at p\=pI. 

We proceed further and define a symmetric function 

n(^Hpl) - \ {^^^HpI - ^e) + P^^\pI + ^e)}. (131) 

It is just this symmetric function on which we will focus in the remainder of the 
section. Let us return to the Minkowski space. This is done by substituting 

p\ — >• exp(— i7r)p^ , u — >• exp(— Z7r)t 

in Eq. ()129p . Furthermore, observing that ImP = Imll, we obtain the following 
representation for the Minkowskian 11: 

^/ X I lmliit)dt „ , , , , 

n(p) = - i : — 4^ + no(p) , (i32) 

where 



7r Jo t — — ie 



2 \pi + p^ — te pi + p^ + te J 

It is necessary to emphasize that e and e are two distinct regularizing parameters, 
e 7^ e. The parameter e is used to regularize the pole at p"^ = —pi, while e sets 
the rules for the branch cut. The most important property of 11 is that the pole 
at = —pi has no imaginary part, by construction. Hence, there is no physical 
particle that corresponds to this pole. In the conventional local field theory the only 
possible additions with no imaginary part are polynomials. Here we encounter a 
new structure which will be discussed in more detail at the end of this section. 

Our goal is to show that a 4D-unitarity-compliant spectral representation holds 
for the Green's function on the brane, at least in the domain where the laws of 4D 
physics are applicable. To this end we turn to the function Q{p,y), defined as 

g = D{p,y)U{p'). (134) 

with the purpose of studying its properties. It is convenient to pass to the momentum 
space with respect to extra coordinates too. Then, the propagator (jl34|) takes the 
form 

H(n2) 

Q{p,q) = , V ■ ■ (135) 
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With these definitions in hand, we can write down the 4D dispersion relation. We 
start from the Kallen-Lehman representation for the propagator (jl35p . As we will 
check below, this representation takes the form 



gfa,),±/^-MMj^^ "°f)-"°W . (136) 

TT JO t — p'^ — It q^^ — — le 



The imaginary part of Q is defined as follows 



1 



Im^(t,g) = -Kdiq^ - t)ReIiit) + \mIiit)V— , (137) 

q^ — t 

where V stands for the principle value of a singular function, 

11/1 1 \ 

'^^^t - 2\q^ -t + i5^ q^ -t-i5)- ^^^^^ 

The fact that Eq. ()13(jj) holds can be checked by substituting ()138j) and ()137j) into 
(|136|) and exploiting the relation 



TT Jo t — p"^ — ie q^ — t 

Ren(g2) - no(g2) + Ii,{p') - Yi{p^) 



(139) 



This turns Eq. (|136p into identity. 

Finally we approach the main point of this section - the dispersion relation for 
G{p,y = 0), the Green's function on the brane. As such, it must have a spectral 
representation with the positive spectral density, as we have already seen from the 
KK-based analysis. The positivity is in one-to-one correspondence with the 4D 
unitarity. 

The dispersion relation can be obtained by integrating (jl36p with respect to q, 



Gip, y = 0) = - r ^ ^^'l^^^ . + Mp') ReDoi-pl) . (140) 

TT JO t — p"^ — le 



According to Eq. ()136|) . the spectral density p is defined as 

pit) = J (141) 

The first term on the right-hand side in Eq. p40|) is conventional while the second is 
not, and we hasten to discuss it. This term has no imaginary part, by construction. 
Hence, it does not contribute to the unitarity cuts in diagrams. Therefore, this term 
does not affect the spectral properties. 
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As was mentioned, in conventional 4D theories only a finite-order polynomial in 

that has no imaginary part can be added to or subtracted from the dispersion 
relation. This is because, normally one deals with Lagrangians which contain only 
a finite number of derivatives, i.e., a finite number of terms with positive powers 
of in the momentum space. In the problem under consideration this is not the 
case, however. In fact, no local 4D Lagrangian exists in our model at all, and yet 
we are studying the spectral properties in terms of the intrinsically 4D variable, 
p^. The theory (jH)) is inherently higher- dimensional because of the infinite volume 
of the extra space. One can try to "squeeze" it in four dimensions at a price of 
having an infinite number of 4D fields. For such a theory there is no guarantee that 
analyticity of the Green's functions in terms of the 4D variable p"^ will hold because 
the effective 4D Lagrangian obtained by "integrating out" the infinite gapless KK 
tower will necessarily contain ,20^ nonlocal terms of the type d~^. (Note that a 
similar prescription for the poles in a pure 4D local theory [HZj is hard to reconcile 
with the path integral formulation [HH]- In our case this is not a concern since the 
theory is not local in four- dimensions in the first place.) 

Therefore, it is only natural that 4D unitarity can be maintained but 4D ana- 
lyticity cannot. Nonanalyticity leads to violation of causality, generally speaking. 
That is to say, the Green's function (jl4(J|) is acausal. Therefore, we have an apparent 
violation of causality in the 4D slice of the entire (4 + A^) dimensional theory which, 
by itself, is causal. The apparent acausal effects can manifest themselves only at the 
scale of the order of ~ 10^^ cm. In fact, as was noted in j2Tj, this is a welcome 
feature for a possible solution of the cosmological constant problem. 

Let us now return to the first term on the right-hand side of Eqs. ()14()j) . Using 
Eqs. fll32|) and ()137|) we can calculate the spectral function which comes out as 
follows: 

Pit) = 2mllmDoit) 

' {tReDo + 2miy + (tlm^o)' ' 



where 



N+2 
2 JV- 



= ^ / (2# ) = (2.rr(N/2) ■ '"3' 

We see that p(t) satisfies the positivity requirement. 

Next we observe that at large momenta, i.e., at p^ Dq{j)) ^ m^, the propagator 
we got has the desired 4D behavior. For the scalar part of the propagator this is 
expected. However, with regards to the tensorial structure this circumstance is less 
trivial. If p^ Do{p) ^ the terms in the braces in Eq. (jl24p . 

T^.T'^'^ - ^TT'I , (144) 

^^For N > 5 the integral in Eq. H14U|) diverges. However, since our model has a manifest UV 
cutoff A, the above integral must be cut off at A. Alternatively, one could use a dispersion relation 
with subtractions. 
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correspond to the exchange of two physical graviton polarizations. Therefore, for 
the observable distances the tensorial structure of the massless 4D graviton (jl44j) is 
recovered. 

On the other hand, for large (super- horizon) distances, Do{p) <^ m^, we get a 
different tensorial structure of the same amplitude, 

jf^.f'^'^ - if f'l . (145) 

This exactly corresponds to the exchange of the six- dimensional graviton. 



4.8 D>6 

Corresponding calculations and results are quite similar to the D = 6 case, with 
minor technical distinctions which we summarize below. For N ^ 2 



I'ab 



2 - TV 
Therefore, we get 



VabK. (146) 



d'h,^ = ^^d^h^. (147) 

Then, the trace of the {/iz^} equation takes the form 

(y) Mil dl K + {N + 2) dl K 

= NT^5'^^\y). (148) 

where = 2 — A^(2 — 36). The above equation which can be used to find the solution 
we are after. We proceed parallel to the six- dimensional case. Let us introduce the 
notation 

KiP.y) = N^g,,{p,y), (149) 

where 

r D{p,y) 

yN = —, 2 r, f \ > — FTaTToT + '^^N horn ■ (loUj 

-k^p^ Dq[p) + u^^ [N + 2) 

The solution of the homogeneous equation takes the form 

e^Nhom = D{p,y)5(-kNp'D^{p) + M^(iV + 2)) . (151) 



Here 
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As in the 6D case, we conclude that that there exists a solution to the equation 



-kNp'^Doip) + u^iN + 2) = Q 

with a complex value of p^. These poles occurs on the nonphysical sheets as long as 
/cat > 0, so the Green's function admits the spectral representation. 

Using the expressions above one readily calculates the tree-level amplitude A, 



u 



N 



^^^""^ - I k^pw, - {2 + N)u^ )r • ^1^^^ 



4.9 6 > (2iV - 2)/3iV 

In this case there are no poles on the physical Riemann sheet. Hence, all the poles 
are of the resonance type. The tensorial structure at large distances is that of the 
D-dimensional theory 

[t.uT'^'' - l^Tf'^ . (153) 

However, the tensorial structure at short distances ^ differs from that of 4D 
massless gravity. Hence, some additional interactions, e.g., repulsion due to a vector 
field, is needed to make this theory consistent with data. 



4.10 b<{2N-2)/3N 

The consideration below is very similar to the 6D case. In perfect parallel with the 
6D case we consider for simplicity only the 6 = case and define the function 

P^n\pI) = ^N^N^2)/{2N-2) -plD,ipE)-^e' ^^^^^ 
which has a spectral representation: 

P<->(pi + «) = i n + . (155) 

TT Jo U - p% pi - Pe - T'f^ 

The residue R is determined by Eq. p30|) while pi is now a solution to the equation 

p"Mp*) = ^ P>0. (156) 

As in the 6D situation, we use the expression p55j) to define a symmetric function 
^^n\pI) - I {p^n\pI - + P^n\pI + ^^)}■ (157) 
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The latter, being continued to the Minkowski space, admits the following spectral 
representation: 



1 lmU{t)dt , (^) 



^n{p) = - / + n'o 'ip) , (158) 

vr JO t — p'^ — le 



where 



nr(i>) - \ ( ,^1 . + ^-4-- ) • (159) 

As previously, e and e are two distinct regularizing parameters, e 7^ e. 
For the Green's function of interest 

D(p,y)n^(p^) 

= 2iV-2 

we repeat the analysis of Sects. 14.51 14.61 and 14.71 to confirm with certainty that the 
function QN{p,y = 0) does admit the spectral representation ()14m) . with a positive 
spectral function, similar to the 6D case, see Eq. (jl42p . 

The expression in Eq. (jl52p interpolates between the four-dimensional and D- 
dimensional patterns. This was already established for the scalar part of the ampli- 
tude. Let us examine the tensorial part. For Do{p) ^ we get 

f^X^'" -]^ff'Y (161) 

This corresponds to two helicities of the 4D massless graviton. In the opposite limit, 
p^ Dq{p) <C m^, we recover the tensorial structure corresponding to the (4 + A^)- 
dimensional massless graviton. 



4.11 No strong coupling problem 

We start from a brief review of the well-known phenomenon - the breakdown of 
perturbation theory for the graviton with the hard mass j36j , occurring at the scale 
lower that the UV cutoff of the theory ^Hl HOI IM] • We then elucidate as to how this 
problem is avoided in the models (jH)), (f73jl . 
The 4D action of a massive graviton is 

Sm = ^J dS^gR{g) + ^ / d^x [h^ - (h^^f) , (162) 

where nig stands for the graviton mass and h^^, = {g^^ — ri^^)/2. The mass term has 
the Pauli-Fierz form jSH] . This is the only Lorentz invariant form of the mass term 
which in the quadratic order in h^^, does not give rise to ghosts [221 • Higher powers 
in h could be arbitrarily added to the mass term since there is no principle, such as 
reparametrization invariance, which could fix these terms. Hence, for definiteness. 
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we assume that the indices in the mass term are raised and lowered by rjniy. Had we 
used g^i, instead, the difference would appear only in cubic and higher orders in h, 
which are not fixed anyway. 

In order to reveal the origin of the problem let us have a closer look at the free 
graviton propagators in the massless and massive theory. For the massless graviton 
we find 



/I 1 1 \ 1 



where 



Vt^u = Vt^u - ■ (164) 

The momentum dependent parts of the tensor structure were chosen in a particular 
gauge convenient for our discussion. On the other hand, there is no gauge freedom 
for the massive gravity presented by the action (jl62j) : hence, the corresponding 
propagator is unambiguously determined, 

X , \ -, (165) 

nig — — le 

where 

f^^,u = - . (166) 

We draw the reader's attention to the 1/m^, 1/m^ singularities of the above propa- 
gator. The fact of their occurrence will be important in what follows. 

It is the difference in the numerical coefficients in front of the rj^yTjap structure 
in the massless vs. massive propagators (1/2 versus 1/3) that leads to the famous 
perturbative discontinuity BHl Ell- matter how small the graviton mass is, 
the predictions are substantially different in the two cases. The structure (jl65p gives 
rise to contradictions with observations. 

However, as was first pointed out in Ref. [22] , this discontinuity could be an 
artifact of relying on the tree- level perturbation theory which, in fact, badly breaks 
down at a higher nonlinear level [SHllini- One should note that the discontinuity 
does no appear on curved backgrounds [HIl E2] ^ another indication of the spurious 
nature of the "mass discontinuity phenomenon." 

To see the failure of the perturbative expansion in the Newton constant Gn one 
could examine the Schwarzschild solution of the model ()162p . as was done in Ref. 
However, probably the easiest way to understand the perturbation theory breakdown 
is through examination of the tree-level trilinear graviton vertex diagram. At the 
nonlinear level we have two extra propagators which could provide a singularity in 
nig up to l/m^. 
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Two leading terms, and I/tti^, do not contribute [lU]; so that the worst 

singularity is l/mg. This is enough to lead to the perturbation theory breakdown. 
For a Schwarzschild source of mass M the breakdown happens 1221 HO] at the scale 

A^~m,(Mm,/M2j)-i/5. 

The result can also be understood in terms of interactions of longitudinal polariza- 
tions of the massive graviton which become strong [S3] . For the gravitational sector 
per se, the corresponding scale reduces to 



rrigimg/Mpi)-^/'' . 

If one uses the freedom associated with possible addition of higher nonlinear terms, 
one can make [S3] the breaking scale as large as 



m,/K/Mpi)i/3_ 

Summarizing, in the diagrammatic language the reason for the precocious break- 
down of perturbation theory can be traced back to the infrared terms in the propa- 
gator (jl65p which scale as 

^ . (167) 
rrig 

These terms do not manifest themselves at the linear level; however, they do con- 
tribute to nonlinear vertices creating problems in the perturbative treatment of 
massive gravity already in a classical theory. 

We will see momentarily that similar problems are totally absent in the propa- 
gator of the model ([H]). For illustrational purposes it is sufficient to treat the N = 2 
case. All necessary calculations were carried out in Sect. 14.41 Therefore, here we 
just assemble relevant answers. 

For iV = 2 and b > 1/3 we find 

2ml - (36- l)j92Do(p) + ie' ^ ' 

In the limit the above expression, as opposed to Eq. (j28j) . is regular. Similar 

calculations can be done in the N > 2 case. The results is proportional to 

P,PuD{p,y) ^^gg^ 



{2 + N)u^ - kNP^Doip) + ie 



which is also regular in the — > limit where it approaches the 4D expression. 
Therefore, we conclude that there is no reason to expect any breaking of perturbation 
theory in the model ^ below the scale of its UV cutoff. 
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If 6 < 1/3 and = 2 we find, by the same token, 



Pt^PuD{p,y) 



2ml - (36- l)p2Do(p) + ie 



(170) 



Again, in tlie limit — > the above expression, in contradistinction with Eq. 
is regular. Moreover, in this limit (and at ?/ = 0) it approaches the 4D expression, 
in a particular gauge. Analogous calculations can be readily done in the N > 2 and 
6 < (2A^ - 2)/3iV case. The results is 



{2 + N)u^ -kNP^ Do{p)+te 



(171) 



This expression is also regular in the rric limit where it arrives at the correct 4D 
limit. We conclude therefore, that in the general case there is no reason to expect 
any breaking of perturbation theory in the model (j75p below the scale of its UV 
cutoff. Note that the expressions (|17Up and (|17ip are singular for small Euclidean 
momenta ~ ""^c- construction this singularity has no imaginary part and 
there is no physical state associated with it. One might expect that this singularities 
will be removed after the loop corrections are taken into account in a full quantum 
theory. These considerations are beyond the scope of the present work. 

An analogy with the Higgs mechanism for non-Abelian gauge fields is in order 
here. For massive non-Abelian gauge fields nonlinear amplitudes violate the unitar- 
ity bounds at the scale set by the gauge field mass. This disaster is cured through 
the introduction of the Higgs field. Likewise, nonlinear amplitudes of the 4D massive 
gravity ()162j) blow up at the scale A^. The unwanted explosion is canceled at the 
expense of introducing an infinite number of the Kaluza-Klein fields in (jSI). 



4.12 Perturbation theory in nonflat backgrounds 

So far we have been discussing perturbations about a brane that has no tension. 
This was done just to learn what kind of gravity was produced by the brane induced 
term on the world-volume. However, the problem of real physical interest is to do 
the same calculation with a brane that has an arbitrary tension. For N = 2 case 
this was studied in Ref. (TOj with the conclusion that the properties obtained in 
Refs. |141 [W\ [7T| lUUj hold unchanged. In the next two subsections we will perform 
the qualitative analysis for > 2 in terms of the Green's function on the brane 
following Refs. jl5[ IHUl , as well as in terms of the KK modes following the method 
of Ref. [71]. We will show that at observable distances the 4D laws of gravity are 
indeed reproduced. 



4.12.1 Propagator analysis 

The nature of gravity on the brane perhaps is simpler understood from the propa- 
gator analysis. The equation for the graviton two-point Green's function (we omit 
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tensorial structures) takes the form 



T5^^\x)M_, (172) 



where T denotes the source (which will be put equal to 1 below) and 

d^+N = -^dAVCG^^dB + higher derivatives, 64 = d^d^,. (173) 
V G 



Using the technique of Ref. jT^I, the scalar part of the solution in the Euchdean 
four-momentum space can be written as 

where D{p,y) is the Euclidean 4-momentum Green's function of the bulk operator 
O^+N, that is 04+ND{p,y) = 6{y)/y'^~^. What is crucial for us is the behavior of 
the Green's function on the brane, 

Let us discuss this expression first. The denominator in p75p consists of two terms. 
The first term, Mpjp^, is what gives rise to 4D behavior. The second term in the 
denominator, M^^^ D^^{p, 0), sets the deviation from the 4D laws and is due to the 
infinite-volume extra bulk. Therefore, in the regime when Mpip'^ dominates over 
M^'^^ D~^{p,0) we get 4D laws, while in the opposite case we obtain the higher- 
dimensional behavior. The question is what is the crossover scale at which this 
transition occurs. To answer this question we need to know the expression for 
D{p,0). Let us start for simplicity with the case when S = 0, i.e., the background 
metric is flat. We will denote the corresponding Green's function by Do{p,y) to 
distinguish it from D{p,y). Moreover, let us drop for a moment higher- derivatives 
in the expression for O^+n- In this case Dq{p, y) is nothing but the Green's function 
of the (4 + A^)-dimensional d'Alambertian. Its behavior at the origin is well known, 

Do{p,y-^0) ~ (176) 

Hence, Doip,Q) diverges and therefore the term M^+^D^^{p,0) = M^+^Do\p,0) 
in Eq. ()175|1 goes to zero. This would indicate that 4D gravity is reproduced on 

^^Note that in the warped case the scale Mi„d differs from that of the flat case by a constant 
multiplier yl^(A). For simplicity this won't be depicted manifestly below. 
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the brane at all distances. However, the UV divergence in p76|) is unphysical. This 
divergence is smoothed out by UV physics [121 EHl ISH Il60j. In reality the bulk 
action and the operator contain an infinite number of high-derivative terms 

that should smooth out singularities in the Green's function in p76|) . Since these 
high-derivative operators (HDO) are suppressed by the scale M*, it is natural that 
the expressions ()176p is softened at the very same scale y ~ . As a result 
one obtains 159* 60^ D^^p^y = 0) ~ Mf-2(1 + 0{p/M^)). Substituting the latter 
expression into (|175j) we find that 

Mr^D-i(p, 0) = M^+'^D.^p, 0) ~ M: . 

Therefore the crossover scale is Vc ~ Mpi/M^ ~ 10^*^ cm. At distances shorter than 
~ 10^^ cm the 4D laws dominate. 

Let us now switch on effects due to a nonvanishing tension S. The background 
in this case is highly distorted. The distortion is especially strong near the brane. 
Let us start again with the case when the HDO's are neglected and contains 
only two derivatives at most. Then, the expression for the D-function reduces to 

D{p,y) ~ Do{p,y)Hya/y), (177) 

where, as before, DQ{p,y) ~ l/y^^^ and JF is some function which is completely 
determined by the background metric (by the functions A, B and C) and JF(0) = 
const. In the region where the solution of the Einstein equations can be trusted, 
can be approximated as follows, J^iyg/y) = {yg/y)°' + c, where a and c are some 
constants determined by A^. If we were to trust this solution all the way down to 
the point y = we would obtain again that M^^^ D~^{p, 0) = and that gravity is 
always four-dimensional on the brane. However, as we discussed above (see also the 
previous section), the existence of high-derivative terms tells us that the background 
solution cannot be trusted for distances y <ti y*. In general, = M~^{ygM^y with 
7 <^ 1 and ^ yg- Thus, for y <^ y^ the higher curvature invariants become large 
in units of M^, and infinite number of them should be taken into account. In order 
to find the effect of this softening, let us take a closer look at the expression (jl77|) . 
There are two sources of singularities in this expression. The first one emerges on the 
right-hand side of ()177p as a multiplier, Dq ~ l/y^~'^; this singularity was discussed 
above in p76|) and is independent of the background geometry. Instead, it emerges 
when the operator 0^+^ is restricted to the quadratic order only. We expect that 
this singularity, as before, is softened at the scale after the higher derivatives 
are introduced in O^+n- Hence, in ()177|) when we take the limit |/ — we should 
make a substitution 

D, ~ Mf-2(1 + 0{p/M..)) . 

On the other hand, the second source of singularity in (jl77|) is due to the function JF. 
This singularity is directly related to the fact that the background solution breaks 
down at distances of the order of y^. As we discussed in the previous section, the 
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UV completion of the theory by HDO's should smooth out this singularity in the 
background solution. In order to get the crossover scale we can use the following 
procedure which overestimates the value of M^+^D~^(p, 0). In the limit y ^ we 
could make the substitution i/g/y Vg/y* in the expression for JF and in ((TTTj). 
Using these arguments we find 



D(p,0)^Mf-^ {yg/y 



+ c 



Moreover, taking into account that ^ yg we get 

Therefore, we conclude that, as in the zero-tension case, the crossover distance^^ in 
the nonzero tension case can be of the order of 10^^ cm. 



4.12.2 Kaluza— Klein mode analysis 

The purpose of this section is to study the effect of a nonvanishing brane tension 
on 4D gravity in terms of the KK modes. The Einstein equations (with up to two 
derivatives) that follow from the action (j^ can be written as 

M2+^ (tZab - IgabTI^ + Ml,6'X6''B6^''\y) [r^, - ]^g^,R 

= Sg,uW'^''\y)- (178) 

Below we consider fluctuations h^y{x, yn) that are relevant to 4D interactions on the 
brane, 

ds^ = A'^{y) ['r]i_,„ + hf,u{x,yn)] dx^dx" 

- B\y)dy' - C\y)y'dnl_,. (179) 

As typically happens in warped backgrounds, equations for graviton fluctuations are 
identical to those for a minimally coupled scalar [721 ES] • The present case is no 
exception. Equation (jl78|) on the background defined in Eq. (|179p takes the form 

M^^a4+iv Ku{x, yn) + O4 KAx, 0) = , (180) 



yN-lA\y) 



where 



= -^OaVGG^'' Ob, 6^ = d^d, . (181) 



^^If we were to assume that the background solution is softened at rather than at y,, we 

would obtain even larger value for the crossover scale. This can be turned around to make the 
following observation. If the background metric softens at M~^, and/or if j/g 3> y*, the value of 
should not necessarily be restricted to 10~^ eV, but can be much higher. Unfortunately, these 
properties does not seem to allow further analytic investigation. 
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To simplify Eq. ()180p we turn to spherical coordinates with respect to ?/„,n = 
1, 2, .., A^, and decompose fluctuations as follows: 

h^,{x,y) = t^,{x)a{y)<p{Sl), (182) 

where the components in Eq. ()182|) satisfy the conditions 



1 Tn?' 
Oi+Ne^,y{x) = —d^d^ef,u{x) = -—e^i,{x), (183) 



Using these expressions we rearrange Eq. ()18U|) as follows: 



' a = —a. (185) 



M^+^y^-^A^{y)j 

Our goal is to rewrite this expression in the form of a Schrodinger equation for 
fluctuations of mass m. We follow the method of Refs. [ZH EH]- It is useful to 
introduce a new function 

X = -^=(^, (186) 

and a new coordinate 

ry B(t) 

u ^ / rfr-ii. (187) 
Jo A[t) 

In terms of these variables Eq. (jl8(J|) takes the form 

\ <f , , A^lil + N-2) m'^Ml.5iy)\ , , , 

where the effective potential V^fr(ti) is defined as 



^eff(«) = -n -TT^ ■ (189) 



du^ \VabJ ■ 

Note that the first two terms in (I188D can be rewritten as 



+ VM -(— + —] (-— + — ^ (190) 
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where 



exp (B) 



(191) 



With appropriate physical boundary conditions the operator on the right-hand side 
of p90|) is selfadjoint positive-semidefinite with a complete set of eigenfunctions of 
nonnegative eigenvalues. 

Let us analyze Eq. ()188p . in particular, the properties of the KK modes following 
from it. What is crucial for our purposes is the value of the KK wave functions on 
the brane, i.e., |x('7i, ?/ = 0)p. The latter determines a potential between two static 
sources on the brane [7T]. We would like to compare the properties oi\x{jn,y = 
which are known jTUinni only for ^ = 0, with the properties obtained at £^ 7^ 0. 

First, we recall the properties of \x{fn^y = 0)p for = 1 and tensionless brane, 
£" = 0. In this case A = B = 1, C = Q and / = 0. Hence, u = y and Kfr(w) = 0. 
Equation (|188j) becomes 



For each KK mode of mass m there is a delta-function attractive potential, the 
strength of which is proportional to the mass of the mode itself. Hence, the higher 
the mass, the more the influence of the potential is. The attractive potential leads 
to a suppression of the wave function at the origin (suppression of \x{jn,y = 0)^). 
Therefore, the larger the mass of a KK state, the more suppressed is its wave function 
at zero. 

Simple calculations in this case yield 



where rc ~ Mpj/M^. This should be contrasted with the expression for |x(m, y = 
0)p in a theory with no brane induced term (i.e., with no potential in Eq. ()192p ). 
In that case |x('^)l/ = 0)^ = 1. We see that the KK modes with masses m ^ 
are suppressed on the brane. The laws of gravity on the brane are provided by light 
modes with m <^ r~^. This warrants ^5 that at distances r <^ rc measured 
along the brane the gravity laws are four-dimensional. 

A similar phenomenon takes place for N > 2, with a tensionless brane. Here 
A = B = C = 1, u = y and Ves{u) = 0. The Schrodinger equation takes the 
form of Eq. (jl88|) with the above substitutions. The total potential consists of an 
attractive potential due to the induced term and a centrifugal repulsive potential. 
Because at y ^ the attractive potential is dominant, one finds properties similar 
to the = 1 case. Heavy KK modes are suppressed on the brane - at distances 
r ^ Tc the brane-world gravity is four-dimensional. The only difference |60j is that 
re ~ Mpi/M^ for A^ > 2. 




(192) 



X(m,y = 0)|2 = 4/(4 + mV2), 
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Let us now turn to the discussion of the case of interest when £^ 7^ and A,B,C 7^ 
1. Here the complete equation (jl88p must be studied. For the solutions that soften 
due to the HDO's close to the brane core we expect that as y — 0, m ~ Hence, 
to study the suppression of the wave functions on the brane one can replace (P /du^ 
in ()188|) by d? /dy"^. The next step is to clarify the role of the potential Kfr(u) that 
is nonzero when we switch on the brane tension £ ^ Since a positive tension 
brane should give rise to an additional attractive potential in space with > 2, we 
expect that ^^(m) is negative at the origin (it should tend to —00 at the origin if 
the HDO's are not taken into account). 

The warp factors A, B and C contain the only dimensionful parameter, Ug. So 
does the potential V^ff(M). Therefore, the maximal value of the potential (if any) in 
the interval < y < i/g should be determined by the very same scale, max{V^fj} ~ 

If the form of the attractive potential were trustable all the way down to small 
values of the coordinate, then an attractive nature of the potential could make easier 
to obtain 4D gravity on a brane as compared to the zero tension case. Unfortunately 
we cannot draw this conclusion since the expression for the potential is not trustable 
below the distance scales y < (see discussions in the previous section) . Although 
y^: is smaller than yg, nevertheless this two scales can have the same order of mag- 
nitude. Based on the discussions in the previous section one should expect that 
the potential in the full theory softens below y^, and does not really give rise to 
a substantial attraction below that scale. On the other hand, the potential could 
give rise to some undesirable results. Indeed, it could produce a bump (a potential 
barrier) at some finite distance from the core somewhere in the interval < y < yg. 
For a parameter range for which this discussion is applicable (i.e., for ^ M*) 
the hight of the bump can be of the order of max{\4ff } ~ . A KK mode with 
the mass m > will not feel the presence of of such Veg. Its wave function will 
have the same properties as in the tensionless brane theory (i.e. the modes with 
m > r'^^ will be suppressed on the brane). However, the wave function of any KK 
mode with the mass m ^ yg^ will be additionally suppressed on the brane because 
of the potential barrier in VeR. The question is whether this effect can alter the 
laws of 4D gravity on the brane at observable distances. If yg is small this effect 
will certainly spoil the emergence of 4D gravity on a brane. The reason is that the 
KK modes that are lighter than y^^ will be additionally suppressed on the brane. 
If these were the "active" modes that participate in the mediation of 4D gravity 
at observable distances in the tensionless case, then having them additionally sup- 
pressed would change the 4D laws. However, if yg is sufficiently large the modes 
which are additionally suppressed are very light m < yg, and, if so, "switching off" 
these modes won't be important for 4D gravity. For instance, if Pg > 10^^ cm, as 
it happens to be the case in the present model, gravity at observable distances will 
not be noticeably different from gravity on a tensionless brane. 

Therefore, we arrive at the following qualitative conclusion. In the worst case, 
gravity on the brane world-volume is mediated by the KK modes that have masses 
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in the band ^ m <^ r^^ . Hence, at distances r <^ i/g the effects of the brane 
tension are negligible and gravity on a brane reproduces the known four- dimensional 
laws. Moreover, in a simple case when ~ r~^, one can think, qualitatively, that 
gravity on the brane is mediated by a 4D graviton of mass 

mg - y;^ - rj^ . 

In the present context this value is of the order of the Hubble scale 

rugr^ H^r^ 10"^^ eV . 
A graviton with such a small mass is consistent with observations. 

5 Brief summary 

We reviewed the arguments why large distance modification of gravity is a promising 
direction toward the solution of CCP. The models that modify gravity at large dis- 
tances emerge in the context of braneworlds with infinite- volume extra dimensions. 
The 4D interaction in these models of brane induced gravity is due to the world- 
volume Einstein-Hilbert term. We discuss in detail the properties of these models 
in dimensions five and higher. The 5D model, although not capable of solving CCP, 
is nevertheless a consistent theory of a large-distance modification of gravity with 
many instructive properties. These properties, that are scattered in the literature, 
were collected in Sect. 3 of the present article. Furthermore, we discussed brane- 
induced gravity in dimension six and higher. These models share certain properties 
with the five-dimensional theory, and, at the same time, they differ from it too. In 
particular, they could lead to solution of CCP. To establish whether or not this is 
the case one needs to perform further detailed calculations the algorithms of which 
were outlined in detail in Sect. 4. 
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